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Preface 


This  research  is  pan  of  the  continuing  development  effort  from  the  Vector  Wave  Equation  (VWE) 
Research  Group.  The  work  reported  in  this  thesis  is  the  algorithmic  manipulation  of  the  VWE  to  form  a 
parallel  computational  approach  which,  when  presented  in  a  graphical  form,  provides  a  specification  tool  for 
Very  Large  Scale  Integration  (VLSI)  implementation. 

It  is  my  firm  belief  that  the  effort  of  the  VWE  Research  Group  can  make  a  significant  contribution 
to  the  "State-of-ihe-Art"  in  the  design  and  analysis  of  structures  that  are  impacted  by  electromagnetic  quan¬ 
tities  or  involve  electromagnetic  metrics  as  a  design  constraint.  As  such,  it  has  been  a  pleasure  and  an 
honor  to  be  involved  in  this  work. 

It  is  appropriate  here,  that  I  acknowledge  and  give  thanks  to  the  people  and  forces  that  made  this 
research  possible.  The  AFIT  computer  environment  is  the  best  that  I  have  ever  been  involved  with.  My 
Thesis  Advisor  Maj  Joe  DeGroat,  encouraged  creative  approaches  to  solutions  rather  than  the  security  of 
classical  methods.  These  factors  combined  to  make  for  a  research  environment  that  was  both  exciting  and 
powerful. 

I  would  be  terribly  remiss  if  did  not  give  special  mention  to  my  wife  Debra,  whose  support, 
encouragement,  and  tolerance  of  the  long  hours  required  to  do  this  work  made  this  thesis  possible. 

Finally,  I  give  thanks  to  the  Holy  Father  and  his  son  Jesus  Christ  in  whom  I  believe. 


Jack  L.  Strauss 
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Abstract 


"An  algorithm  has  been  specified  for  the  hardware  implementation  of  the  numerical  solution  of  the 


electromagnetic  fields  of  an  arbitrary  current  source.  The  algorithm,  defined  as  the  Vector  Wave  Equation 

,  \  >  C  '\ 

(VWE),  solves  for  the  magnetic  (ff)  and  electric  fields,  as  well  as  the  vector  potentiaI*(A)  of  a  finite 


length  arbitrary  current  source.  The  specified  algorithm  has  been  verified  through  FORTRAN  simulation 


to  produce  results  accurate  to  within  2  decimal  places  for  a  500  sub-element  dipole. 


The  VWE  forms  a  model  which  is  algorithmically  symmetric  with  respect  to  the  cartesian  coordi¬ 


nate  system.  As  such,  the  VWE  lends  itself  to  highly  parallel  and  concurrent  computational  techniques. 


This  property  of  the  algorithm  makes  it  an  excellent  candidate  for  implementation  by  a  Very  High  Speed 


Integrated  Circuit  (VHSIC)  class  processor.  Investigation  of  a  parallel,  highly  concurrent  architectural 


implementation  has  yielded  preliminary  results  that  computational  savings  of  a  factor  of  3  and  a  throughput 


rate  increase  of  5  orders  of  magnitude  is  attainable. 


This  research  has  shown  that  an  application  specific  VHSIC-class  processor  array  has  sufficient 


computing  power  to  support  interactive  calculation  of  a  set  of  equations  which  solve  for  the  magnetic  and 


electric  fields,  as  well  as  the  vector  potential  of  an  arbitrarily  current  source. 
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ARCHITECTURAL  IMPLICATIONS  OF 


A  PARALLEL  COMPUTATIONAL  APPROACH 
TO  THE  VECTOR  WAVE  EQUATION 

I.  Introduction 

This  thesis  is  the  study,  development,  and  analysis  of  a  parallel  computational  approach  to 
the  set  of  electromagnetic  field  equations  defined  below  as  the  Vector  Wave  Equations  (VWE).  The  charter 
of  this  effort  was  to  specify  a  parallel  processing  engine  for  the  VWE  that  could  be  implemented  in  either 
software  or  hardware.  As  such,  the  derivation  of  a  parallel  algorithm,  specification  of  computational  flow, 
and  an  analysis  of  the  algorithm's  behavior  with  respect  to  error/accuracy  considerations  became  the  focus  of 
this  work. 

This  study  was  motivated  by  the  the  end  objective  of  developing  an  interactive  CAD/CAM  envi  - 
ronment  for  aerodynamic  systems  design  where  a  structure's  electromagnetic  observability,  specifically  its 
radar  cross  section,  is  known  early  in  the  design  phase  of  a  project.  The  capability  of  having  a  structure's 
observability  known  and  interactively  updated  as  the  design  progresses  and  is  changed  requires  calculation 
of  near-  field  electromagnetic  wave  equations.  Numerical  solutions  to  these  equations  is  a  computationally 
intensive  problem  that  requires  millions  of  floating  point  calculations. 

At  present,  software  approaches  on  classical  computer  architectures  have  extremely  long  run  times. 
Because  of  this  constraint,  only  simplified  antenna  configurations  and  geometries  are  examined.  Conse  - 
quently,  there  is  a  need  for  the  development  of  a  very  fast  customized  VHSIC-class  processor  engine  for 
calculating  the  radiated  fields  of  arbitrary  antennas.  This  processor  engine  would  be  part  of  a  workstation 
which  integrates  the  processor's  electromagnetic  computing  capabilities  along  with  interactive  aerodynamic 
structural  design  techniques,  such  that  electromagnetic  metrics  become  an  integral  part  of  the  design  phase. 

The  design  phase  is  the  most  effective  time  to  influence  the  Radar  Cross  Section  (RCS)  of  an 
object.  Once  the  design  has  been  implemented,  RCS  reduction  must  rely  on  the  addition  of  heavy  radar 
absorbing  material,  and  active/passive  cancellation.  Radar  absorbing  material  helps  to  dissipate  the  incident 
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energy  upon  its  surface,  thereby  reducing  the  RCS  of  the  target  However,  this  material  can  add  con  - 
siderable  weight  and  thereby  degrade  the  operational  characteristics  of  the  structure.  Active/passive  can  - 
cellation  techniques  that  reduce  the  RCS  becomes  intractable  after  design  implementation  due  to  the  large 
size  of  the  structures  and  current  analysis  capabilities.  Consequently,  it  is  more  cost  effective  to  make 
structural  changes  within  a  CAD  environment  during  the  design  phase. 

With  the  limitations  to  simplified  geometries  for  numerical  analysis,  complex  structures  such  as 
an  aerospace  vehicle  require  RCS  measurements  with  scale  models  at  compact  ranges,  or  with  full  size 
models  at  outdoor  ranges.  This  approach  is  time  consuming,  costly,  and  subject  to  scaling  problems  and 
poor  weather  conditions.  However,  the  incorporation  of  a  VHSIC -class  processor  engine,  capable  of  high 
speed  calculation  of  the  electromagnetic  equations  into  an  interactive  CAD  workstation,  makes  possible  real 
time  observations  of  electromagnetic  metrics.  As  such,  the  impact  of  modifications  to  the  object's  shape 
and  composition  material's  on  its  RCS  are  known  during  its  design. 

This  introductory  text  will  include  a  discussion  of  the  electromagnetic  background  and  development 
of  the  VWE.  The  approach  used  in  the  derivation  of  the  parallel  algorithm  will  be  stated.  Motivation 
and  basic  assumptions  in  the  derivation  of  the  error  expression  will  be  discussed  and  referenced.  Finally, 
some  thoughts  on  the  utilization  and  usefulness  of  the  results  of  this  study  will  be  touched  on. 

Electromagnetic  Background 

The  basis  for  modern  antenna  theory  lies  with  the  set  of  formulas  known  as  Maxwell’s  equations. 
By  utilizing  Maxwell's  equations,  the  total  radiation  pattern  can  be  derived  for  a  given  antenna  (acting  as  a 
scatterer).  The  total  field  of  an  antenna  is  the  sum  of  the  incident  and  the  scattered  fields  of  that  antenna. 

To  find  the  scattered  field,  the  incident  field  is  subtracted  from  the  total  field.  Once  the  current  distribution 
on  the  scatterer  is  found,  the  radiation  integral  given  by  Maxwell  can  be  numerically  computed. 

A  mid-point  summation  technique  was  employed  to  arrive  at  a  numerical  solution  for  the  vector 


potential  (A)  [Jones,  1985:6].  Initially  it  was  anticipated  that  the  vector  potential  would  be  used  as  an 
intermediate  result  in  obtaining  the  final  electric  (E)  and  magnetic  (H)  fields.  The  solution  for  the  vector 
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potential  involves  an  approximation  where  the  dipole  is  divided  into  a  finite  number  of  sub-elements. 


Equation  (1)  defines  the  mathematical  formula  relating  the  vector  potential  solution  to  the  general  antenna 


problem  as  depicted  in  Figure  1.  The  evaluation  of  the  complex  integral  in  equation  (1)  is  necessary  to 


analyze  even  the  simplest  antenna.  The  desired  magnetic  and  electric  fields  are  derived  from  equations  (2) 


and  equation  (3).  The  mathematical  procedure  is  a  vector  curl  operation,  for  which  this  study  uses  cartesian 


coordinates. 


where 


where 
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However,  it  was  found  that  the  numerical  result  for  the  vector  potential  is  not  necessary  in 
computing  the  magnetic  and  electric  fields.  The  mid-point  summation  method  employed  to  obtain  the 
vector  potential  can  be  used  for  the  magnetic  and  electric  fields  directly  [Hoyt,1986].  Equations  (4)  and  (5) 
below  show  the  magnetic  and  electric  fields  in  the  "x”  direction  for  cartesian  coordinates.  Equations  for  the 
y  and  z  vector  components  are  of  equal  complexity  and  similar  in  form  to  those  in  the  x  direction  for  the 
general  antenna  configuration.  Evaluation  of  these  equations  is  further  complicated  by  current  densities  (J) 
being  complex  quantities  in  the  general  case. 


where  h  =  length  of  the  dipole  per  number  of  sub-elements 


r  =  I  r  -  rj  I  =  sqrt  {(x-xi)2+(y-yi)2+(z-zi)2). 

Hoyt  derived  and  validated  these  equations,  defined  here  to  be  the  discrete  Vector  Wave  Equations 
(VWE),  with  two  current  densities  of  different  geometries  on  a  simple  z-directcd  dipole.  The  algorithm  was 
evaluated  with  varying  numbers  of  sub-elements  to  determine  the  number  required  to  give  an  accuracy  of 
two  decimal  places.  The  result  of  Hoyt's  study  show  that  the  numerical  summation  algorithm  produces 
results  accurate  to  within  1%  for  95  percent  of  the  observation  points  for  a  500  sub-element  dipole.  This 
study  exploits  the  inherent  symmetry  of  the  algorithm  with  respect  to  the  cartesian  coordinate  system  and. 
as  such,  leads  to  a  highly  parallel  and  concurrent  computational  approach. 
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Approach 

With  Hoyt's  results  verifying  that  an  acceptable  level  of  accuracy  can  be  obtained  with  500  sub- 
elements,  and  the  knowledge  that  any  meaningful  construct  will  be  formed  from  thousands,  possibly 
millions  of  finite  length  dipoles,  clearly  the  critical  path  to  reducing  the  computational  intensity 
is  to  calculate  the  quantities  within  the  summation  formula  as  efficiently  as  possible.  Therefore 
this  thesis  effort  is  concerned  with  the  i^1  calculation. 

The  approach  used  in  the  derivation  of  the  computational  expressions  which  yield  the  desired 
results  is  threefold.  First,  the  technique  of  a  top  down  functional  decomposition  was  performed  for  each 
of  the  vector  quan-tities  of  interest.  This  procedure,  yielded  a  database  of  computational  expres¬ 
sions  in  which  intermediate  results  could  be  studied  throughout  the  calculation  process.  Second,  the 
calculation  flow  of  each  vector  quantity  was  mapped  graphically  to  form  the  construct  of  a  data¬ 
flow  graph.  These  data-flow  graphs  yield  visual  evidence  of  the  inherent  symmetry  of  these  expres¬ 
sions  and  provide  the  parallel  computational  specification  desired  as  a  result  of  this  work.  Third,  the 
data-flow  graphs  were  overlayed  to  find  the  levels  of  concurrency  that  are  attainable. 

The  processes  of  functional  decomposition  and  the  mapping  into  the  directed  graphs  is  detailed 
in  chapter  two  of  this  thesis.  Results  in  computational  requirements  reduction  at  various  levels  of 
computational  concurrency  is  detailed  in  chapter  four,  the  analysis  section  of  the  thesis.  Finally,  a 
complete  set  of  calculation  flow  graphs  are  presented  in  appendix  B  for  reference  and  concise  access. 

Error  Expression 

There  are  two  types  of  error  associated  with  any  computer  computation.  The  first  type  is  error 
attributable  to  errors  in  data.  This  type  of  error  is  a  function  of  the  inherent  error  in  finite  word  length  and 
data  representation.  Also,  the  error  transmitted  as  the  value  of  some  function  to  another  process  creates  a 
cumulative  input  error  to  subsequent  computations.  The  second  type,  generated  error,  is  a  grouping  of 
factors  such  as  truncation,  algorithmic  cancellation  of  significant  digits,  and  others  that  are  generated  by  the 
particular  implementation  of  the  function  [Cody  and  Waite,  1980:  11].  Of  the  two  types,  the  programmer 


4 


v 


v 

* 

■  i 

t1 


.V 


or  designer  has  direct  control  of  generated  error  but  must  observe  the  systems  performance  to  gain  insight 
into  the  "effects"  and  "behavior"  of  inherent  and  transmitted  error. 


i 

The  analysis  of  the  cumulative  effects  of  transmitted  error  requires  extensive  simulation  for  / 

I 

nontrivial  systems.  Analytic  techniques  rely  on  the  Stochastic  Model  [Taylor,  1983:  391 1  where  an 
arbitrary  input  is  processed,  and  the  result  is  of  the  form  Y[X(n)]  plus  some  delta  representing  the  error 
indicative  of  the  process.  Manually  tracking  these  types  of  data  for  a  system  of  even  moderate  size  is  error 
prone,  at  best.  The  problem  becomes  intractable  for  complex  systems. 

Chapter  three  of  this  thesis  implements  the  concept  of  the  Stochastic  Model  in  terms  of  inherent 

I 

and  generated  error  to  gain  visibility  of  the  behavior  of  the  Parallel  VWE  with  respect  to  cumulative  j 

transmitted  error.  The  process  is  performed  in  MACSYMA  and  is  presented  at  a  tutorial  level. 

Utilization  of  Results 

In  this  final  section  of  introductory  text,  a  brief  discussion  of  the  intended  utilization  and  resultant 
usefulness  of  the  results  of  this  thesis  effort  is  appropriate.  Three  areas  will  be  addressed:  the  target 
processing  environment,  intended  use  of  the  Error  Expression,  and  the  use  of  MACSYMA  in  handling 
large  sets  of  equations. 

Target  Processing  Environment.  As  stated  by  way  of  the  end  objective  being  a  VHSIC-class 
processor  engine  for  incorporation  into  a  CAD  work-station,  the  initial  target  processing  environment  is  a 
hardware  implementation.  However,  the  results  of  this  study  have  yielded,  by  means  of  the  data-flow 
graphs,  several  levels  of  processing  environment  independent  computation  specifications  that  can  and  should 
be  implemented  and  tested  in  software.  Exploration  into  alternate  architecture  software  implemcntauon  of 
any  of  the  levels  of  computational  concurrency  is  sure  to  yield  time  savings  far  greater  than  any  classical 
computation  techniques. 

Error  Expression.  The  intended  use  of  the  derived  error  expression  discussed  above  is  as  an 
analytical  model  for  hardware  implementation  design  trade-off  decisions.  Because  the  level  of 
abstraction  is  purely  independent  of  any  elementary  function  implementation,  the  result  is  a  model 


that  a  VLSI  designer  can  use  without  constraint  to  any  one  particular  elementary  function 
implementation.  This  type  of  model,  when  complemented  with  empirical  data  from  simulation  or 
known  levels  of  accuracy  from  previous  designs,  can  aid  in  the  assurance  of  overall  system  accuracy 
being  known  during  the  design  phase. 

Use  of  MACSYMA.  One  of  the  most  daunting  aspects  of  the  analysis  of  nontrivial  sets  of 
equations  is  maintaining  analytic  clarity  and  accuracy  of  the  expression  as  the  set  of  equations  is 
manipulated  from  one  step  to  the  next  MACSYMA,  a  symbolic  equation  processor,  has  proven 
through  the  course  of  this  study  to  be  a  tool  capable  of  handling  the  enormous  book-keeping  tasks  as 
well  as  providing  a  full  range  of  mathematical  functions  to  handle  any  computational  task  required.  The 
use  of  MACSYMA  in  this  theses  effort  is  extensive  and  fully  documented.  The  intent  of  the  presen- 
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The  process  of  functional  decomposition  of  the  VWE  equations  through  the  use  of  MACSYMA  is 
presented  in  a  step  by  step  explanation  of  the  Parallel  Vector  Potential  derivation.  The  results  of  the  more 
complex  expressions  for  the  electric  and  magnetic  fields  are  given  in  a  less  tutorial  form  for  brevity.  The 
mapping  of  the  parallel  algorithm  into  the  data-flow  graph  construct  is  also  shown  below. 

The  process  is  one  of  substituting  variables  into  an  equation.  These  variables  represent 
"allowable"  elementary  functions.  Allowable  is  connoted  to  mean  either  a  hardware  functional  element  or  a 
software  process  implemented  in  the  target  processing  environment. 


Figure  2  Decomposition  Algorithm 


The  variables  arc  coded  such  that  computational  and  algorithmic  clarity  is  maintained.  That  is,  the 
coding  scheme  should  clearly  define  computational  levels  in  the  overall  calculation,  while  maintaining  end- 
to-end  algorithmic  correctness.  For  this  study  the  ALPHA  characters  A,  B,  C, ...  denote  distinct 
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computational  levels  while  the  concatenated  numeric  designators  1, 2, 3, ...  uniquely  define  a  specific 
elementary  function.  For  example,  ”E5"  might  mean  a  specific  multiplication  at  the  "E"  computational 
level.  This  process  is  achieved  by  the  algorithm  given  above  in  Figure  2. 

Vector  Potential  Derivation 

The  cartesian  expression  for  the  Discrete  Vector  Potential  is : 


=  ^  (A  X  +  Ay  + 


x  +  «y  + 


where 


Ax  =  G  *  (  jxri  +  jxii  ) 

Ay  =  G  *  (  jyri  +  jyii  ) 

Az  =  G  *  (  jzri  +  jzii  ) 


noting  that  jxri  is  the  real  part  of  the  complex  current  density  in  the  x 
direction,  jxii  is  the  imaginary  part  with  similar  notation  being  used  for  the 
other  coordinate  system  elements, 


which  is  the  free  space  form  of 
Greens  Function, 


r  is  the  vector  distance  between  the  ith  element  and  the 
observation  point;  P  is  the  wave  number  (27t/X). 

These  expressions  are  entered  into  MACSYMA  as  follows: 

g:(%EA(-%l*b*r))/(4*%pi*r); 

ax:g*(jxri+%l*jxii); 

ay:g*(jyri+%l*jyii); 

az:g*(jzri+%l*jzii); 


The  resulting  MACSYMA  presentation  for  the  three  vector  coordinate  system  elements  arc  shown 


below: 


-  %i  b  r 


(jxri  +  %i  jxii)  %e 


4  %pi  r 
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-  %i  b  r 


(jyri  +  %i  jyii)  %e 

ay;  . 

4  %pi  r 

-  %i  b  r 

(jzri  +  %i  jzii)  %e 

az;  - 

4  %pi  r 

Once  entered,  the  top-down  functional  decomposition  of  the  x  directed  vector  potential  may  begin. 

The  process  is  identical  for  the  y  and  z  coordinate  system  elements. 

From  the  above  expression  for  Ax  (ax)  the  real  and  imaginary  parts  are  extracted.  The  real  and 

imaginary  parts  of  the  Vector  Potential  in  the  x  direction  are  defined  as  axr  and  axi  below: 

realpart(ax) 

axr;  jxii  sin(b  r)  +  jxri  cos(b  r) 

4%pi  r 

Imagpart(ax) 

axi;  jxii  cos(b  r)  -  jxri  sin(b  r) 

4  %pi  r 

These  real  and  imaginary  parts  may  be  operated  on  according  to  the  above  substitution  algorithm. 
First,  however,  recall  that  "r"  is  the  vector  distance  between  the  i^1  element  and  the  observation  point.  As 
the  discrete  VWE  given  here  in  cartesian  coordinates,  the  cartesian  expression  for  r  is  required.  This  subsu  - 
tution  gives  an  initial,  fully  decomposed  form  of  the  real  and  imaginary  parts  (denoted  by  axrO  and  axiO)  as: 

axrO;  2  2  2 

(jxii  sin(b  sqrt((z  -  zi)  +  (y  -  yi)  +  (x  -  xi)  )) 

2  2  2 

+  jxri  cos(b  sqrt((z  -  zi)  +(y-yO  +(x-xi)  ))) 

2  2  2 
/(4  %pi  sqrt((z  -  zi)  +(y-yi)  +(x-xi)  )) 
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2  2  2 

axiO;  (jxii  cos(b  sqrt((z  -  zi)  +(y-yi)  +(x-xi)  )) 

2  2  2 

-  jxri  sin(b  sqrt((z  -  zi)  +(y-yi)  +(x-xi)  ))) 

2  2  2 

/(4  %pi  sqrt((z  -  zi)  +  (y  -  yi)  +(x-xi)  )). 

Following  the  development  process,  the  First  substitution  is  Al,  A2,  and  A3  for  (x-xi),  (y-yi),  and 
(z-zi).  At  this  point  in  the  substitution  process,  these  are  the  only  allowable  elementary  operations  capable 
of  being  performed.  All  other  calculations  within  the  equation  are  contingent  upon  the  completion  of  these. 
Therefore,  there  is  a  delineation  of  a  computational  level.  The  status  of  the  equation  after  this  first  level  of 
substitution  is: 


We  note  two  things  here.  First,  the  quantity  4jr  (4  %pi)  is  considered  as  an  input  constant  such 
that  a  separate  multiplication  is  not  required.  Secondly,  in  the  third  level  of  computation  shown  below,  the 
function  Cl  is  a  triple  summation.  The  development  of  the  VLSI  Triple  Summer  will  be  a  future  thesis 
project  and  was  accepted  as  an  allowable  elementary  function  for  this  study.  With  the  triple  summer,  the 
status  after  the  third  level  of  substitution  is: 

jxii  sin(b  sqrt(Cl))  +  jxri  cos(b  sqrt(Cl )) 

axr3;  - 

4%pi  sqrt(CI) 


and  axi3; 


jxii  cos(b  sqrt(CI))  -  jxri  sin(b  sqrt(Cl)) 


4%pi  sqrt(CI) 


The  square  root  operation  is  performed  next  via  substitution  of  D1 .  The  status  after  the  forth  level 


jxii  sin(b  01)  +  jxri  cos(b  D1) 


4  %pi  D1 


and  axi4; 


jxii  cos(b  D1 )  -  jxri  sin(b  D1 ) 
4  %pi  D1 


Next,  the  multiplication  of  (5  and  D1  is  performed,  for  which  we  substitute  El .  Here  we  also  have 
the  multiplication  of  the  input  constant  An  (4  %pi)  and  D1  replaced  with  E2.  The  status  after  the  fifth  sub  - 
stitulion  is: 


jxiisin(EI)  +  jxricos(EI) 


and  axi5; 


jxii  cos(EI)  -  jxri  sin(EI) 


/ 


The  process  continues  with  the  substitutions  FI,  F2,  and  F3  for  cos(El),  sin(El),  and  the  inver  - 
sion  process  of  E2,  respectively.  Note  that  an  inversion  process  is  used  rather  than  a  division  to  provide  a 
more  general  specification  for  the  hardware  implementation.  The  status  after  level  six  is: 

axi6;  (jxii  F2  +  jxri  FI)  F3 

and  axi6;  (jxii  FI  -  jxri  F2)  F3  . 

The  next  substitutions  are  for  the  multiplications  of  FI  and  F2  with  the  complex  components  of 
the  current  density.  The  status  after  the  seventh  level  of  computation  is: 

axr7;  F3(G4  +  G1) 

and  axi7;  F3(G3-G2). 

The  next  operations  to  be  performed  are  the  remaining  addition  and  subtraction  from  above.  The 
status  after  level  eight  is: 
axr8;  F3  HI 

and  axB;  F3  H2  . 

Finally,  we  multiply  the  quantities  above.  Recall  that  F3  is  the  resultant  inversion  quantity  from 
level  six.  The  final  results  are: 

axr9;  J1  The  real  part  of  the  vector  potential  in  the  "X"  direction 

axi9;  J2  The  imaginary  part  of  the  vector  potential  in  the  “X“  direction. 

Mechanically  the  process  is  one  of  scanning  the  the  equation  to  check  if  an  allowable  elementary 
function  can  be  performed.  If  one  can,  it  is  coded  by  the  conventions  stated  above.  That  is,  the  function  is 
given  the  ALPHA  designator  for  the  current  computational  level  followed  by  a  numeric  designator  uniquely 
defining  that  specific  function.  The  resultant  coded  variable  is  substituted  into  the  equation,  and  the  same 
process  is  continued  until  no  more  elementary  functions  can  be  performed.  The  ALPHA  designator  is 
incremented,  and  the  process  is  continued  until  a  single  variable  remains  signifying  that  the  calculation 


process  is  complete. 


Note  that  for  this  study,  the  above  process  was  accomplished  manually.  However,  clearly  this 
process  is  one  that  lends  itself  to  and  would  benefit  by  automated  techniques.  Also  the  letter  "I"  was  not 
used  in  defining  computational  levels  to  avoid  confusion  between  alpha/numeric  characters. 

The  results  of  this  process  provide  us  with  a  database  of  computational  expressions  at  each  level 
where  the  substitution  designators,  Al,  B3,  etc.,  indicate  floating  point  operations  required  in  the  evaluauon 
of  the  equation.  Once  this  process  is  accomplished  for  each  of  the  coordinate  system  elements,  a  listing  of 
the  required  substitutions  is  generated.  The  listing  for  the  entire  cartesian  set  is  generated  by  eliminating 
the  redundant  calculations.  The  combined  calculations  form  a  parallel  computational  approach  to  the  vector 
potential.  Below  is  a  listing  for  the  discrete  calculation  of  the  real  part  of  Ax  and  the  combined  cartesian 


Level  nine: 


18! 

! 
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subst("H2"1"G3"-"G2'',%); 
subst(’,H3"1,'G5”+"G8",%); 
subst("H4","G7"-"G6",%); 
subst("H5","G9”+"Gl2",%); 
subst("H6","Gl  1  "-"G10",%); 


ft 


subst("J1H,HF3 . H1M,%);  subst("J1".“F3"*"H1",%); 

subst("J2","F3"*"H2",%); 

subst(,,J3","F3H*"H3",%); 

subst("J4","F3"*"H4'',%); 

subst("J5”rF3"*"H5",%); 

subst("J6","F3"*"H6",%); 


v, 


«* 


s 


It  is  noted  that  there  is  a  one  to  one  correspondence  of  the  substitution  designators  to  the  nodes  of 
a  data-flow  graph.  A  graphical  representation  is  generated  by  mapping  the  above  listings  into  data-flow 
graphs.  The  data-flow  graphs  for  the  real  part  of  the  "x"  directed  vector  potential  and  the  combined  cartesian 
set  are  shown  in  Figures  3  and  4. 


.v 

The  substitution  process  for  the  vector  potential  is  presented  in  MACSYMA  syntax  in  Appendix 


A.  Data-flow  graphs  for  all  vector  quantities  are  given  in  Appendix  B. 


* 
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Magnetic  Field  Derivation 


The  equation  for  the  Discrete  Magnetic  Field  in  cartesian  coordinates  is: 

H  =  X  (Hx  +  Hy  +  Hz) 

i 

The  complexity  of  the  equations  precludes  a  tutorial  presentation.  The  processing  of  these  field 
equations  is  of  such  magnitude  that  manipulation  becomes  intractable  without  an  equation  processor  such  as 
MACSYMA.  For  example,  at  the  point  when  the  cartesian  expression  for  the  vector  distance  ”r"  is 
substituted  into  the  expression  for  the  real  part  of  the  magnetic  field  in  the  "x"  direction  forming  "HxrO", 
the  equation  has  the  form: 

2  2  2  2 
((cos(b  sqrt((z  -  zi)  +(y-yi)  +  (x-xi)  ))  (jyri  (z  -  zi)  -  jzri  (y  -  yi))  +  sin(b  sqrt((z  -  zi) 

2  2  2  2  2 
+  (y  -  yi)  +  (x  -  xi)  ))  (jyii  (z  -  zi)  -  jzii  (y  -  yi)))  /sqrt((z  -  zi)  +  (y  -  yi)  +  (x  -  xi)  ) 

2  2  2 

-  b  (cos(b  sqrt((z  -  zi)  +  (y  -  yi)  +  (x  -  xi)  ))  (jyii  (z  -  zi)  -  jzii  (y  -  yi))  -  sin(b 
2  2  2 

sqrt((z  -  zi)  +  (y  -  yi)  +  (x  -  xi)  ))  (jyri  (z  -  zi)  -  jzri  (y  -  yi)))) 

2  2  2 
/(4  %pi  ((z  -  zi)  +  (y  -  yi)  +(x-xi)  )). 

The  substitution  process  for  the  magnetic  field  is  carried  out  similarly  to  that  of  the  vector  poten  - 
ual  described  above.  The  listings  are  generated  and  the  data-flow  graphs  are  formed.  The  data-flow  graphs 
lor  the  real  part  of  the  "x"  directed  magnetic  field  and  the  combined  cartesian  set  is  given  below  in  Figures  5 
and  6.  The  complete  substitution  process  for  the  magnetic  field  is  presented  in  MACSYMA  syntax  in 
Appendix  A.  Data-flow  graphs  for  all  vector  quantities  are  given  in  Appendix  B. 


Electric  Field  Derivation 

The  equation  for  the  Discrete  Electric  Field  in  cartesian  coordinates  is: 

E  =  (Ex  +  Ey  +  E^j 
i 

The  complexity  of  the  electric  Held  equations  is  even  greater  than  the  equations  for  the  magnetic 
field.  For  example,  at  the  same  point  illustrated  above,  when  the  cartesian  expression  for  the  vector  dis  - 
tance  "r"  is  substituted  into  the  expression  for  the  real  part  of  the  electric  field  in  the  "x"  direction  forming 
"ExiO",  the  equation  has  the  form: 


2  2  2 
(sin(b  sqrt((z  -  zi)  +(y-yi)  +(x-xi)  ))  { 


2  1 

+  (jxii  (b - ) 

2  2  2 
(2  -  zi)  +(y-yi)  +  (x-xi) 


3  b  (jzri  (x  -  xi)  (z  -  zi)  +  jyri  (x  -  xi)  (y  -  yi)) 


2  2  2 
sqr1((z  -  zi)  +  (y  -  yi)  +  (x  -  xi)  ) 


3bfxri  2  2 

- )  ((z  -  zi)  +  (y  -  yi)  ) 

2  2  2 
sqrt{(z  •  zi)  +  (y  -  yi)  +  (x  -  xi)  ) 
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3  2 

+  ( - b  )  (jzii  (x  -  xi)  (z  -  zi)  +  jyii  (x  -  xi)  (y  -  yi))) 

2  2  2 
(z-zi)  +  (y  -  yi)  +(x-xO 


2  2  2  3b(jzii(x-xi)(z-zi)+jyii(x-xi)(y-yi)) 

+  cos(b  sqrt((z  -  zi)  +  (y  -  yi)  +  (x  -  xi)  ))  ( - 

2  2  2 
sqrt((z  -  zi)  +  (y  -  yi)  +  (x  -  xi)  ) 


2  2  2 
sqrt((z-zi)  +  (y-yi)  +(x-xi)  ) 


2 

+  jxri(b  -  — 


1  2  2 

- ))((z-zi)  +  (y  -  yi)  ; 

2  2  2 
(z-zi)  +  (y  -  yi)  +(x-xi) 


3  2 

+  ( - b  )  (jzri  (x  -  xi)  (z  -  zi)  +  jyri  (x  -  xi)  (y  -  yi)))) 

2  2  2 
(z-zi)  +  (y  -  yi)  +(x-xi) 


2  2  2  32 

/(4%pi((z-zi)  +(y-yi)  +{x-xi)  )  ) 


The  substitution  process  for  the  electric  field  is  carried  out  similar  to  that  for  the  vector  potential 
and  magnetic  field  described  above.  The  listings  are  generated,  and  the  data-flow  graphs  are  formed.  The 
data-flow  graphs  for  the  real  part  of  the  "x"  directed  electric  field  and  the  combined  cartesian  set  is  given 
below  in  Figures  7  and  8.  The  complete  substitution  process  for  the  electric  field  is  presented  in 
MACSYMA  syntax  in  Appendix  A.  Data-flow  graphs  for  all  vector  quantities  arc  given  in  Appendix  B. 
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It  is  important  to  note  that  the  data-flow  graphs  developed  here  provide  several  types  of 
specification  data.  A  programmer  can  use  these  graphs  as  the  algorithmic  specification  to  reduce  processing 
by  means  of  eliminating  redundant  computations.  The  computer  architect  will  note  the  high  degree  of 
concurrent  processing  that  can  be  achieved  in  an  "Alternate  Architecture"  approach.  Further,  the  VLSI 
designer  now  has  a  system  architecture  specification  at  the  functional  element  level.  The  impact  and 
implications  of  this  parallel  computational  approach  is  developed  in  chapter  four  of  this  thesis. 


In  this  section  we  seek  to  develop  an  analytic  expression  which  will  provide  a  model  for  imple  - 
mentation  trade-off  analysis  for  the  parallel  VWE  with  respect  to  accuracy/error  considerations.  In  doing  so 
we  define  an  algorithmic  entity  as  the  VWE  Core  Equation.  This  expression  is  common  to  all  three  vector 
quantities  of  interest  to  this  thesis  effort  The  development  of  the  error  expression  is  a  bottom-up  process 
in  which  error  terms  are  assigned  to  each  of  the  allowable  elementary  functions  and  to  the  input  data.  The 
resultant  terms  are  substituted  into  the  VWE  core  and  the  computation,  as  specified  in  chapter  two,  is  car  - 
ried  oul  The  final  expression  includes  the  arithmetic  quantity  desired  and  an  upper  bound  on  the  cumulauve 
results  of  the  functional  and  generated  error.  This  resultant  expression  forms  a  model  for  analysis  of  the 
VWE  Core  Equation  with  error  terms  as  the  parametric  variables. 

The  text  below  is  organized  into  the  definition  of  the  VWE  core  and  assignment  of  the  error  terms 
for  the  allowable  elementary  functions,  detailed  presentation  of  the  bottom-up  derivation  of  the  error  ex¬ 
pression,  and  the  resulting  error  model  which  can  be  used  for  implementation  trade-off  analysis. 

The  VWE  Core 

We  see  from  the  data-flow  graphs  that  the  VWE  Core  Equation  forms  two  quanuucs: 

Cos(  P  sqrt(  (x-xi)^  +  (y-yi)2  +  (z-z.i)^  ) 
and 

Sin(  P  sqrt(  (x-xi)^  ♦  (y-yi)"  +  (z-zi)^  ) 

For  this  discussion  we  will  consider  the  cosine  term  only.  Graphical  representation  of  the  VWE  Core 
Equation  is  shown  below  in  figure  9. 


s  \ 


Data  Representation 

eO 

Addition  /Subtr  action 

•  1 

Multiplication 

e2 

Square  Root 

e3 

Sin /Cos 

e4 

Inversion 

eS 

Triple  Summation 

e6 

Cumulative  Error  (Level :  L) 

egL 

Table  1  ErrorTerms 


These  error  terms  are  defined  as  the  total  error  inherent  in  the  input  data  ( bit/word  length),  and  the 
error  generated  as  the  function  value  (resolution,  truncation,  etc.).  A  notational  convenience  will  be  used 
when  the  cumulative  expression  detracts  from  analytical  clarity.  It  is  given  a  designator  defining  its  comp  - 
utational  level  and  is  treated  as  a  separate  arithmetic  quantity  at  higher  levels.  This  approach  facilitates  im  - 
plementation  of  the  stochastic  error  model  for  complex  systems  and  yields  a  final  error  expression  that  is 
independent  of  any  specific  elementary  function  implemen-tation.  Thus,  the  error  expression  provides  only 
specification  data  and  does  not  restrict  the  VLSI  designer  to  any  one  (analyzed)  implementation. 

Bottom-uo  Derivation  of  Error  Model 

As  in  chapter  two,  this  section  will  be  tutorial  with  respect  to  the  level  of  MACSYMA  syntax 


used  in  its  presentation.  The  intent  is  to  provide  an  example  to  future  users  for  similar  applications. 


% 

k 

i 

! 

We  start  by  entering  the  core  into  MACSYMA: 
core:cos(b*sqrt(  (x-xi)A2  +  (y-yi)A2  +  (z-zi)A2 )); 

(c24)  core; 

2  2  2 

(d24)  cos(b  sqrt((z  -  zi)  +(y-yi)  +(x-xi)  )) 

We  proceed  similarly  to  the  approach  taken  in  the  calculation  flow  analysis.  The  text  below  il  - 
lustrates  the  steps  taken  in  this  analysis.  The  following  series  of  substitutions  replace  the  input  values  x, 
xi,  y,  ...zi  with  a  construct  of  the  form  x  +  eO,  xi  +  eO, ...  zi  +e0.  The  "eO"  term  accounting  for  the  in  - 
herent  error  of  the  finite  length  binary  representation  used  in  any  computer  architecture. 

subst(x-xi+2*e0.x-xi,%) ; 

2  2  2 
oos(b  sqrt((z  -  zi)  +  (y-yi)  +(-xi  +  x  +  2e0)  )) 

subst(y-yi+2*e0,y-yi,%); 

2  2  2 
cos(b sqrt((z  -  zi)  +(-yi  +  y  +  2e0)  +(-xi  +  x  +  2e0)  )) 

subst(z-zi+2*e0.z-zi,%) ; 

2  2  2 
cos(b  sqrt((-  zi  +  z  +  2  eO)  +  (-  yi  +  y  +  2  eO)  +  (-  xi  +  x  +  2  eO)  )) 

Next  we  perform  the  arithmetic  operations  of  level  "A"  and  add  the  associated  error  terms  for  the 

subtraction  operation  "el”.  We  note  here  the  procedure  that  will  be  used  in  the  derivation  of  the  error  ex  - 

pression.  To  form  an  upper  bound  expression,  we  assume  the  form  of  the  general  addition  of  error  terms. 

That  is,  we  will  not  attempt  to  arithmetically  cancel  error  terms  at  the  atomic  level.  To  do  so  would  be  to 

assume  more  about  a  target  architecture  than  is  intended  by  this  study.  As  such,  the  following 

MACSYMA  substitution  processes  are  of  the  form  (x-xj)  =>  (x-Xj)  +  el. 

substfAl”  +  el+2*e0,x-xi+2*e0,%)$ 
subst("A2"  +  e1+2*e0,y-yi+2*e0,%)$ 
subst("A3"  +  el  +2*e0,z-zi+2‘e0,%); 
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The  resulting  form  of  the  VWE  Core  after  the  first  level  of  comp-utation  is  significantly  more 
complex  than  that  of  chapter  two,  again  exemplifying  the  need  to  use  an  equation  processor.  The  status  is 
shown  as  corel  below: 
corel: 

2  2  2 
cos(b  sqrt((A3  +  el  +  2  eO)  +  (A2  +  el  +  2  eO)  +  (AI  +  el  +  2  eO)  )) 

Next  we  multiply  the  inner  terms  to  enable  use  to  proceed  We  can  also  see  the  effect  of  cum  - 
ulative  error,  noting  the  "3er  and  "12e0"  terms. 
expand(%); 

2  2  2 
cos(b  sqrt(A3  +  2  el  A3  +  4  eO  A3  +  A2  +  2  el  A2  +  4  eO  A2  +  AI 

2  2 
+  2  el  Ai  +4eOAl  +3  el  +  I2e0el+I2e0  )) 

The  calculations  of  level  two  are  performed  by  first,  substituting  for  the  operation.  Then,  wc 
substitute  in  the  error  terms  associated  with  the  squaring  operation  (same  as  multiplication),  the  process  is 
shown  below: 

subst(',Br,"Ai"A2t%)$ 

Subst("B2"1”A2”A2,%)$ 

subst(”B3","A3"A2,%); 

C0S(b  sqrt(B3  +  B2  +  B1  +  2  el  A3  +  4  eO  A3  +  2  el  A2  +  4  eO  A2  +  2  el  Ai 

2  2 

+  4eOA1+3e1  +I2e0el+12e0  )) 

subs»("Br  +  e2."Bi",%)$ 
substfB2’  +  e2,"B2",%)$ 
subst("B3"  +  e2,"B3",%); 

The  status  of  the  computation  at  level  two  is: 

core2; 

cos(b  sqrt(B3  +  B2  +  B1  +  2  el  A3  +  4  eO  A3  +  2  el  A2  +  4  eO  A2  +  2  el  Ai 

2  2 

+  4eOAl +3e2  +  3e1  +I2e0el  +  12e0  )) 
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To  perform  the  computation  of  level  three  (the  triple  summation),  we  need  to  substitute  the  entire 
inner  expression  as  MACSYMA  will  not  allow  sub-atomic  substitutions.  The  process  and  the  status  at 
level  three  is  given  below: 

subst("C1"+e6+2*e1  +"A3"+4*eO*"A3"+2*e1  *"A2”+4*eO*“A2"+2*e1*"A1  "+4*eO*"A1"+3*e2+3* 
el  A2+12*eO*e1+12*eOA2,’B3"  +  "82”  + 

"81  “+2*e1  *"A3"+4*eO**A3’+2*e1*"A2"+4*eO*"A2"+2*er"A1  "+4*eO*“A1  "+3*e2+3*e1  A2+1 2*e 
0*e1+12*e0A2,%); 

Note,  we  have  added  the  error  term  e6  as  part  of  the  substitution. 

core3; 

cos(b  sqrt(Cl  +  66  +  4e0A3  +  A3  +  2e1A2  +  4e0A2  +  2e1A1  +  4  eO  A1 

2  2 

+  3 e2  +  3 el  +12e0e1 +2e1 +12e0  )) 

At  this  point  we  wish  to  simplify  the  expression.  We  will  define  "eg3”  as  the  generated  error  at 
level  three.  The  development  is  as  follows: 

let  X  *  +  e6  +  4  eo  A3  +  A3  +  2  el  A2  +  4  eO  A2  +  2  el  A1  +  4  eO  A1 

2  2 

+  3e2  +  3e1  +  I2e0e1 +2e1  +  12e0 

this  expression  represents  all  of  the  terms  inside  the  square  root  operator  except  Cl  Recall  the  identity  : 

sqrt(A+B)  =  sqrt[(A+B)*(A/A)j  -  sqrt  [(A+B)/A)  *  sqrt  (A) 

we  then  define  eg3  =  sqrt[(Cl+X)/Cl]  and  perform  the  substitution  with  resulting  expression  shown 
below: 

subst(b*sqrt("Cl")*eg3,  b*sqrt("Cl“  +e6+2*eU"A3"+4*e0*"A3"+2*el*"A2’ 
+4‘eO*“A2"+2*6l*"A1"+4*eO*"Ar+3*e2+3‘e1  A2+12*eO*e1  +  l2*eOA2),%). 

core3.l ; 

cos(b  eg3  sqrt(Cl)) 
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We  noie  here  that  eg3  has  a  multiplicative  relationship  to  the  desired  atomic  expression  "sqrt(C  1 )" 
As  we  proceed  to  level  four,  we  have  not  violated  any  arithmetic  axiom  as  eg3  has  the  form  of  the  square 
root  operation  inherent  in  it  as  shown  in  the  development  above. 

The  computation  and  status  at  level  four  is  given  below: 

subst(b‘("DV  +  e3  )*  eg3,b*eg3*sqrt("Cl")  ,%); 
core4; 

cos(b  eg3  (D1  +  e3)) 

We  can  expand  the  inner  expression  as  follows: 

core4  1 ; 

cos(b  eg3  Dl  +  b  e3  eg3) 

At  this  step  we  see  that  the  expression  yields  a  multiplicative  error  term  which  is  not  readily 
separable  from  within  the  argument  of  the  cosine  term.  As  such,  we  proceed  to  level  Five  with  the 
cumulative  error  as  shown  above.  The  operation  is  performed  and  the  error  term  is  substituted  below. 


II 


V. 

V 


A 


subst(”E1".b*"D1",%); 

subst(“Er+e2,"E1",%): 


This  last  substitution  yields: 


cos[  El  eg3  +  (  p  e3  eg3  )  +  e2  ] 

which  by  trigonometric  substitution  has  the  form: 

cos(El  eg3)  cos((  (3  e3  eg3  )  +  e2)  -  sin(El  eg3)  sin((  p  e3  eg3  )  +  e2) 

We  sec  that  at  this  stage  that  we  can  only  get  the  desired  term: 

COS(EI) 


if  and  only  if: 


eg3  «=  1 

and 


(  p  e3  eg3  )  +  e2  =  n2n  ,  n=  0,1, 2, 3. 
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We  now  have  a  closed  form  expression  that  in  a  limited  sense  models  the  behavior  of  the  VWE 
calculation.  With  the  exception  of  the  inversion  or  division  function  we  now  have  an  expression  which 
includes  all  of  the  arithmetic  operations  (and  the  associated  error  terms)  required  in  the  calculation  of  the 
VWE. 

The  VWE  Error  Model 

The  VWE  Error  Model  follows  from  the  derivation  above.  We  saw  that  the  desired  value  "cos(El )" 
is  attained  if  and  only  if: 

eg3  «  1 
and 

(P  e3  eg3  )  +  e2  =  r\2n  ,  n=  0,1, 2, 3,... 

where  eg3  =  sqrt[(Cl+X)/Cl]  =  sqrt[1  +  X/C1] 

Cl  «{x-xi)2  +  (y-yi)2  +  (z-zi)2 

*  (x2  -  2  x  xi  +  xi2  +  y2  -  2  y  yi  +  yi2  +  z2  -  2  z  zi  +  zi2) 

X  -  4  eO  A3  +  A3  +  2  el  A2  +  4  eO  A2  +  2  el  A1  +  4  eO  A1  +  e6 

2  2 
+  3e2  +  3el  + 12  eO  el  +2  el  +  12  eO 

A1  »  x-xi 
A2  =  y-yi 
A3  =  z-zi 

These  relationships  constitute  a  model.  The  designer  with  known  error  metrics  for  the  funcuonal 
elements  which  implement  the  required  elementary  operations  can  substitute  those  metrics  into  the  expres  - 
.>ions  and  gain  visibility  of  the  impact  that  any  particular  implementation  will  have  on  the  system's  accu  - 
racy  well  before  bread  boarding  a  design.  This  process  supports  design  trade -o I  f  decisions  with  respect  to 


"system  accuracy"  specifications  at  a  much  earlier  and  cost  effective  phase  of  the  project. 


It  is  imperative  to  note  that  this  model,  alone,  may  not  be  suitable  for  such  design  decisions. 
Exhaustive  simulation  as  well  as  derivation  of  a  family  of  VWE  error  expressions  and  models  may  be 
required.  The  complexity  of  the  VWE  and  the  many  possible  architectural  implementations  will  drive  this 
issue. 

The  model  developed  in  this  study  is  a  baseline.  Further  development  and  analysis  is  beyond  the 
scope  of  this  thesis. 


IV.  Analysis 


In  this  section,  the  total  computational  requirements  of  the  parallel  VWE  is  analyzed  and  the 
calculated  savings  at  various  levels  of  computational  concurrency  is  presented.  The  data-flow  graphs 
developed  in  chapter  two  are  used  to  extract  significant  data  for  this  analysis.  The  significant  data  derived 
from  the  graphical  representation  includes  visibility  of  operational  redundancy  within  the  coordinate 
components  of  a  specific  vector  quantity  and  across  the  three  vector  quantities  of  interest.  A,  E,  and  H. 
Visibility  of  data  dependencies  and  maximal  concurrency  possible  in  the  calculation  of  all  components  of 
A,  E,  and  H  is  readily  available.  Finally,  the  graphs  provide  a  means  for  comparison  of  implementation  of 
the  equation  at  various  levels  of  computational  concurrency. 

Removal  of  Operational  Redundancy 

Classically,  the  real  and  imaginary  parts  of  each  of  the  vector  components  would  be  discretely 

calculated  without  regard  to  the  fact  that  several  of  the  intermediate  results  are  reusable  in  the  calculation  of 
other  coordinate  system  elements.  For  example,  when  calculating  the  real  part  of  Hx  ,  there  are  30  floating 

point  operations  involved.  There  are,  noting  symmetry,  30  floating  point  operations  involved  in  the  calc  - 
ulation  of  the  imaginary  part  of  Hx  .  However,  as  clearly  visible  in  the  graphical  representation,  26  of 

these  operation  are  identical.  Therefore  both  quantities  may  be  obtained  in  34  flops.  Further  analysis  of  the 
graphical  representation  shows  that  the  complete  vector  component  set  for  both  the  real  and  imaginary  parts 
of  the  magnetic  field,  H,  can  be  calculated  in  just  74  flops  as  apposed  to  180  when  computed  classically. 
The  analysis  is  summarized  in  table  2  for  all  vector  quantities  of  interest 

Also  shown  in  table  2  is  an  entry  for  "E  &  H  fields”.  This  entry  represents  the  removal  of  all 
computational  redundancy  within  the  complete  E  and  H  graphical  representations  and  is  the  overlaying  of 
the  2  computational  processes  to  form  one  computational  process.  The  savings  of  26  flops  achieved  by 
this  overlaying  is  due  to  a  common  core  of  operations  present  in  these  vector  quantities  similar  to  that 
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Table  2  Computational  Requirements 


defined  in  chapter  three  as  the  VWE  Core  Equation.  Because  of  this  common  core  of  operations,  the  vector 
potential.  A,  can  be  integrated  into  the  fully  parallel  calculation  of  the  three  vector  quantities  of  interest, 
such  that  all  results  are  obtainable  in  202  flops  as  apposed  to  594  flops  when  computed  classically, 
representing  nearly  66%  computational  savings.  The  impact  of  these  results  are  shown  in  graphical  form  in 
figure  10  where  the  computational  requirements  are  plotted  for  a  given  number  of  dipole  elements. 
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Figure  10  Number  FLOPS  vs  Dipole  Elements 
for  a  given  Computational  Approach 


Brora 


Maximal  Concurrency 


Concurrency  in  this  thesis  is  understood  in  an  architectural  sense.  That  is,  concurrency  rr  cans  the 
simultaneous  execution  of  a  number  of  discrete  operations.  Specifically,  in  this  thesis,  concurrency  means 
the  simultaneous  execution  of  floating  point  arithmetic  operations. 

The  visibility  into  the  data  dependencies  that  the  data-flow  graphs  provides  allows  us  to  observe  12 
distinct  computational  levels  when  the  entire  set  of  vector  quantities  are  overlayed  as  depicted  on  the  fol¬ 
lowing  page  in  figure  1 1. 


This  represents  a  98%  reduction  in  linear  time  when  compared  to  the  classical  approach.  However 
this  dramatic  time  reduction  is  attainable  only  if  the  processing  environment  supports  concurrent  floating 
point  arithmetic  operations  at  the  maximum  width  of  the  VWE  calculation.  The  maximum  operational 
width  of  the  VWE  is  36,  occurring  at  the  seventh  computational  level,  G.  These  36  operations  include  30 
floating  point  multiplies  and  six  adder/subtractors. 

Architectural  Considerations 

With  the  end  objective  of  developing  an  interactive  CAD/CAM  environment  for  aerodynamic 
systems  design,  the  target  processing  environment  is  assumed  to  be  a  VHSIC  class  processor,  that  is,  a 
hardware  implementation.  However,  the  computational  specification  data  available  from  the  parallel 
approach  is  target  processing  environment  independent  in  that  any  current  or  future  architecture  can  realize 
processing  gain  by  implementing  the  parallel  VWE  at  the  maximum  level  of  computational  concurrency 
supportable  by  that  architecture. 

As  shown  in  the  table  2  above,  discrete  calculation  of  the  i111  vector  element  requires  594  FLOPS. 
Assuming  that  the  accuracy  verified  by  Hoyt  is  sufficient,  we  will  have  500  of  these  1th  calculations  per 
dipole  element.  If  we  define  the  execution  time  of  the  average  floating  point  arithmetic  operation  to  be  Ta, 

then  we  derive  an  expression  to  approximate  the  time  required  for  the  solution  of  the  discrete  VWE  which 
we  define  as  T{jjscrete: 

"’"discrete  =  [500  *  594Ta]  *  N, 

where  N  =  the  number  of  dipole  elements  in  the  structure. 

Using  the  results  of  the  parallel  approach  where  all  redundant  compulations  have  been  eliminated, 
as  shown  in  table  2,  there  are  202  FLOPS  required.  Defining  Tcombmed  to  be  time  required  lor  this 

solution  we  have: 


"’"combined  =  (500  *  202Ta  ]  *  N. 


This  approach  when  compared  to  the  discrete  approach  provides  a  speed  up  of  nearly  3.  It  is 
significant  to  note  that  this  speedup  can  be  achieved  on  the  same  architecture  as  that  of  the  discrete  by 
simply  recoding  the  algorithm  to  use  all  partial  results.  The  speedup  actually  achieved  on  a  given 
architecture  may  be  somewhat  less  than  3  due  to  the  width  of  the  computation  and  the  resulting  need  to 
store  partial  results  in  memory  versus  registers. 

If  however  the  the  given  architecture  supports  concurrent  floating  point  arithmetic  then  the 
parallelism  shown  in  the  graphical  representation  may  be  exploited.  In  the  fully  concurrent  algorithm  there 

are  12  distinct  computational  levels.  We  define  the  average  execution  time  of  the  concurrent  floating  point 
arithmetic  process  to  be  Tg  con  .  The  resultant  expression  for  the  time  required  to  calculate  the  parallel 

VWE,  defined  to  be  T  ...  is: 

parallel 


T  „  ,  =  [500*  12T  „  ]  *  N. 

parallel  1  a-corr 


Comparing  with  T^^ 


gives: 


T ..  .  T 

discrete  a 

* 

- =  49  5  - 

Tparaflel  Ta-con 


this  speed  up  of  49.5  can  be  achieved  if  the  architecture  supports  concurrent  floating  point  operations  at  the 

maximum  width  of  the  parallel  VWE  and  the  ratio  T  to  T  is  =  1 .  The  maximum  width  is  36.  If  less 
v  a  aeon 

than  36  functional  units  are  available,  the  speedup  would  be  proportionally  reduced. 

A  block  diagram  for  a  possible  architecture  that  is  defined  as  the  hardware  implementation  of  the 
fully  overlayed  data-flow  graph  of  figure  1 1  is  shown  below  in  figure  12.  The  architecture  depicted  in 
Figure  12  is  one  in  which  there  is  a  functional  unit  for  each  node  in  the  data-flow  graph.  The  functional 
units  arc  hardwired  so  that  as  results  become  available  they  are  passed  to  the  next  unit.  Due  to  the 
computation's  data  dependencies,  processor  functional  units  at  all  but  the  currently  active  level  sit  idle. 

This  results  in  very  low  component  utilization. 
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Figure  12  Block  Diagram  and  Utilization 
of 202  Funtional  Unit  Architecture 


Further  improvements  in  speedup  can  be  achieved  by  applying  pipelining  and  other  architectural 
techniques.  The  limit  of  these  approaches  would  be  a  hardwired  pipelined  architecture  capable  of  producing 
the  complete  vector  component  set  every  major  cycle.  A  block  diagram  of  this  architecture  is  given  in 
figure  13  below. 

The  depicted  pipelined  architecture  requires  the  same  202  functional  units  as  the  previous 
architecture  with  the  addition  of  the  interstage  registers  and  pipeline  control  hardware.  The  improvement  in 
system  performance  is  due  to  a  result  being  available  every  major  cycle  (once  the  pipe  is  full)  as  menuoned 


earlier  and  a  significant  increase  in  component  utilization  when  a  large  number  of  dipole  elements  arc  to  he 

v" 


34 


r. 


calculated.  Utilization  figures  for  this  architecture  are  not  static  as  before  but  are  dependent  on  the  number 
of  dipole  elements  to  be  input  into  the  pipe.  The  limit  of  component  utilization  approaches  100%  for  large 
numbers  of  dipoles  to  be  calculated  and  the  curve  is  shown  in  figure  13. 


The  above  architecture's  performance  can  be  expressed  as  follows: 

Tpiplined“(500  *  Tq^.)  *  N, 


where  Tmc  is  the  time  of  a  major  cycle.  This  expression  shows  that  the  speed-up  is  limited  to  approx  - 


imately  two  orders  of  magnitude  for  a  single  processor  system. 


A  means  of  achieving  further  speed-up  is  through  the  application  of  multiple  processors  each 
capable  of  computing  the  VWE  component  set  and  taking  advantage  of  the  architectural  considerations 


presented  thus  far. 

An  architecture  capable  of  supporting  30  floating  point  multiplies,  18  floating  point 
additions/subtractions,  a  square  root,  sin/cos  and  an  inverse  or  division  operation  may  compute  the  complete 
ith  VWE  calculation  in  12  computational  levels.  This  is  an  architecture  with  51  functional  units  and  a 
block  diagram  is  given  in  figure  14.  These  51  functional  units  are  reconfigurable  to  implement  the  various 
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Figure  14  Block  Diagram  and  Utilization 
of  51  Funtional  Unit  Architecture 


computational  levels  specified  in  figure  11  and  hardwired  in  the  previous  architectures.  This  architecture 
will  require  extensive  use  of  registers  and  control  circuitry  that  will  increase  complexity  and  area  con  - 
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siderations.  However,  because  of  the  significantly  reduced  number  of  functional  units,  component  util  - 
ization  factors  are  improved  as  shown  in  figure  14.  Implementation  of  this  architecture  utilizing  VHSIC 
technology  and  possible  wafer  scale  integration,  yields  a  processing  engine  that  may  be  incorporated  as  the 
processor  of  a  multiprocessing  array  system. 

The  capabilities  of  the  above  system  is  shown  in  the  definition  of  Tconcurrent  below: 

^concurrent  “  [500  *  2Tvhsid  *  N 

_  i 

Number  of  Processors 

where  Tvj,sjc  is  the  execution  time  of  the  average  arithmetic  operation  on  the  VHSIC  processor  engine. 

Comparing  an  array  of  the  above  processors  to  the  classical  single  processor  architecture  mentioned 
previously  results  in  a  speedup  (S)  of: 

Ta 

S  =  49.5  * _ *  Number  of  Processors  . 

Tvhsic 

We  now  seek  an  expression  that  bounds  S  in  terms  of  the  number  of  processors  used  in  the  array. 

A  lower  bound  given  by  assuming  one  processor  and  that  the  ratio  T^T^^c  =  1.  That  is,  a  VHSIC 

implementation  of  the  functional  units  are  of  the  same  order  of  speed  as  the  average  execution  time  of  the 
classic  architecture.  The  resultant  expression  is: 

Ta 

49.5  5  S  <  49.5  * _ *  Number  of  Processors  . 

Tvtisc 

The  upper  bound  makes  no  assumption  as  to  Ta/Tv[jSjc  or  number  of  processors.  Considering  a 
conservative  lower  limit  of  2  for  Ta/Ty^^  and  rounding  the  scaler  on  the  right  hand  side  of  the  inequality 

to  50  we  can  express  the  equation  as: 

49.5  <  S  <  100  *  Number  of  Processors  . 

To  achieve  a  desired  5  orders  of  magnitude  improvement,  the  above  expression  indicates  that  an 
array  of  1000  VHSIC  processor  engines  are  required. 


A  realistic  evaluation  of  an  aerodynamic  structure  may  easily  require  1  million  dipole  elements. 
Using  the  discrete  approach  on  a  Cray  1,  this  computation  would  require  8.25  hours  assuming  that  the  100 
MFLOP  rate  could  be  uniformly  sustained.  Even  if  the  parallel  algorithm's  results  are  utilized  to  remove 
all  redundant  operations,  the  computation  would  still  require  2.8  hours.  This  result  clearly  does  not  support 
an  interactive  environment.  However,  the  computation  for  these  same  1  million  dipole  elements  would 
require  only  0.6  seconds  on  the  1000  processor  array  as  described  above.  This  is  a  processing  rate  capable 
of  supporting  an  interactive  environment 

The  performance  of  the  various  architectures  discussed  above  are  compared  with  respect  to  time  in 
figure  15  below.  For  this  comparison,  the  specific  architecture  time  variables  are  all  assumed  to  be  10-ns. 
The  effect  of  this  generalization  is  nominal  with  respect  to  the  magnitudes  calculated  below. 


Architectural 

Approach 


NUMBER  DIPOLE  ELEMENTS 

♦  Classical  ^discrete 

0  Classical  w/redundency  removed 
^combined 

■  Parallel  Processor  Tpara|lej 
o  Pipelined  Processor  Tpjpe)ined 

*  Array  (1000)  Tconcurrent 


Figure  15  Time  vs  Dipole  Elements 
For  a  given  Architectural  Approach 
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This  study  was  motivated  by  the  the  end  objective  of  developing  a  CAD/CAM  environment  for 
aerodynamic  systems  design.  In  particular,  the  capability  of  having  a  structure’s  electromagnetic  observ  - 
ability  known  and  interactivly  updated  as  the  design  progresses.  Current  analytical  approaches  require 
calculation  of  near-field  electromagnetic  wave  equations.  Numerical  solutions  to  these  equations  is  a 
computationally  intensive  problem  that  requires  millions  of  floating  point  calculations. 

This  thesis  represents  the  study,  development,  and  analysis  of  a  parallel  computational  ap  - 
preach  to  the  set  of  electromagnetic  field  equations  defined  as  the  Vector  Wave  Equations  (VWE).  This 
effort  has  produced  an  initial  specification  for  a  parallel  processing  engine  for  the  VWE  that  can  be  imple  - 
mented  in  either  software  or  hardware.  The  analysis  points  excitingly  to  a  conclusion  that  an  interactive 
environment  with  respect  to  electromagnetic  metrics  is  attainable.  The  cost  of  this  capability  is  a  complex 
array  of  VHSIC  class  processors.  Further  cost  must  be  considered  to  address  the  following  issues. 

The  perspective  of  this  development  effort  was  that  a  model  of  the  structure  existed  in  a  form  com  - 
patible  with  the  cartesian  representation  of  the  VWE.  As  such,  the  issue  of  creating  that  model  must  be  ad  - 
dressed.  That  task  will  likely  be  of  a  similar  magnitude  as  this  forward  processing  task  has  been.  The  issue 
of  the  model  with  respect  to  it's  shear  size  must  be  grappled  with.  Finally,  the  grim  task  of  I/O  handling 
and  post  processing  must  be  completed  in  order  to  attain  a  truly  interactive  design  environment 

Other,  possibly  more  attractive,  issues  are  those  of  alternate  applications  for  the  VWE.  An 
example  might  be  the  study  of  it's  utility  in  long  lead,  high  frequency  layout  analysis  where  the  wave 
length  of  the  frequency  is  a  significant  fraction  of  the  conduction  path.  Further,  the  applicability  of  these 
results  to  a  general  class  of  related  partial  differential  equations  may  make  finite  element  calculations 
practical  for  heat  flow,  fluid  dynamics,  and  other  "vector  intensive’’  areas  of  study. 

Final  recommendations  for  follow-up  research  with  respect  to  this  thesis  are  in  the  areas  of 
empirical  verification  of  the  computational  time  savings  and  hardware  implementation.  It  was  stressed 
throughout  this  text  that,  although  the  objective  was  a  hardware  implementation,  the  algorithmic 


specification  of  the  data-flow  graphs  should  be  coded  and  tested  in  software.  This  area  would  include 
development  of  a  VWE  simulator  at  various  levels  of  computational  concurrency  and  testing  the  algorithms 
performance  on  systems  such  as  the  Cray.  Further  software  projects  should  include  a  simulator  of  the  VWE 
with  the  allowable  elementary  functions  parameterized  by  function  resolution  and  word  size.  This  project 
will  give  empirical  visibility  into  the  algorithms  performance  with  respect  to  accuracy  considerations. 

Hardware  implementation  issues  include  time-area  analysis,  optimal  scheduling  and  utilization 
research,  and  VHSIC  technology  concerns.  The  form  of  the  Parallel  VWE,  as  developed  here,  may  lead  to 
implementation  by  wafer  scale  integration.  Further  investigation  into  alternate  architectures  combined  with 
the  resolution  of  some  of  the  above  issues  may  likely  be  the  path  to  attaining  a  very  "state-of-the-art"  goal 
of  an  interactive  design  environment 


This  appendix  provides  the  supplementary  material  tor  the  develop¬ 


ment  of  the  Parallel  VWE.  Presented  here  is  the  full  development  process  of 


the  three  vector  quantities  (A,  H,  and  E)  in  MACSYMA  syntax  The  text  is 


organized  into  three  sections.  The  sections  are  the  development  of  the 


Parallel  Vector  Potential,  the  Parallel  Magnetic  Field,  and  the  Parallel  Electric 


Field.  The  individual  sections  are  further  divided  into  the  presentation  of  the 


computational  flow  which  is  the  development  of  the  parallel  electromagnetic 


quantity  and  the  command  listing  which  implements  that  process  These 


command  listings  are  what  is  mapped  into  the  data-flow  graphs  The  data¬ 


flow  graphs  are  presented  in  Appendix  B. 


This  section  illustrates  the  calculation  analysis  of  the  Vector  Potential 


The  following  pages  will  show  the  step  by  step  procedure  of  the  computation 


of  the  Parallel  Vector  Potential 


The  vector  potential  in  the  "x"  direction  will  be  the  example  through 


level  four  with  the  complete  complex  six-tuple  shown  for  the  remainder  of 


the  process  The  text  is  in  MACSYMA  syntax  and  presented  as  the  user  might 


see  it  on  the  screen 


( c 2 2 )  ax. 


-  751  b  r 


(jxn  +  Si  jxii)  fie 


A  Spi  r 


realnarttax)  The  real  part  is  extracted  with  MACSYMA. 


(c23)  axr; 

(d23) 


jxii  sm(b  r)  +  jxri  cos(b  r) 


4  Spi  r 


imaqpart(ax)  The  imaginary  part  is  extracted  with  MACSYMA 
(c24)  axi, 


(d24) 


Jxii  cos(b  r)  -  jxri  sin(b  r) 


4  Spi  r 


Here  we  beoln  the  computation  process 

The  example  below  recalls  the  real  part  of  Axr,  then  substitutes  the  expanded  version  of  "r"  into  the 
expression. 

AxrO;  Substitute  the  expanded  version  of  ”r"  into  the  expression. 

(c25)  axrO. 

2  2  2 

(d25)  (jxii  sin(b  sqrt((z  -  zi)  *  (y  -  yi)  +  (x  -  xi)  )) 

2  2  2 

+  jxri  cos(b  sqrt((z  -  2i)  +  (y  -  yi)  +  (x  -  xi)  ))) 

2  2  2 

/(4  3pi  sqrt((z  -  zi)  +  (y  -  yi)  +  (x  -  xi)  )) 


AxlO;  Substitute  the  expanded  version  of  “r"  into  the  expression. 
(c26)  axiO, 

2  2  2 

(d26)  (jxii  cos(b  sqrt((z  -  zi)  +  (y  -  yi)  (x  -  xi)  )) 

2  2  2 

-  jxri  sin(b  sqrt((z  -  zi)  +  (y  -  yi)  +  (x  -  xi)  ))) 

2  2  2 

/( 4  ftp!  sqrt((z  -  zi)  +  (y  -  yi)  (x  -  xi)  )) 


After  the  first  level  of  computation: 


Axrl; 

(c27)  axrl, 


Axil; 


(c23)  axil, 


2  2  2  2  2  2 
jxii  sin(b  sqrt(A3  +  A2  +  A1  ))+  jxri  cos(b  sqrt(A3  +  A2  +  A1  )) 


2  2  2 

4  Spi  sqrt(A3  +  A2  +  A1  ) 


2  2  2  2  2  2 
jxii  cos(b  sqrt(A3  +  A2  +  A1  ))  -  jxri  sm(b  sqrt(A3  +  A2  +  A1  )) 


2  2  2 

4  Xpi  5qrt(A3  +  A2  +  A 1  ) 


After  the  second  level  of  computation: 


Axr2; 


(c35)  axr2; 


jxii  sin(b  sqrt(B3  +  B2  ♦  BU)  +  jxri  cos(b  sqrt(B3  +  B2  ♦  BD) 
4  Spi  sqrt(B3  ♦  B2  ♦  B 1 ) 


Axi2; 


(c47)  axi2; 


jxii  cos(b  sqrt(B3  +  B2  +  B I ))  -  Jxri  stn(b  sqrt(B3  +  B2  +  B 1)) 


4  Xpi  sqrt(B3  +  B2  +  B 1 ) 


After  the  third  level  of  computation: 

Axr3; 

(c38)  axr3; 


(d38) 


jxii  sin(b  sqrt(C  1 ))  +  jxri  cos(b  sqrUC  I )) 


4  Spi  sqrtfC  1 ) 


Axi3; 

(c39)  axi3, 

(d39) 


jxil  cos(b  sqrt(C  1 ))  -  jxri  sin(b  sqrt(C  1 )) 


4  Kpi  sqrt(CI) 


After  the  forth  level  of  computation: 


After  the  seventh  level  of  computation: 

Axr7; 

(c203)  axr7; 

(d203)  F3 (64  +  61 ) 

Axi7; 

(c204)  axi7; 

(d204)  F3 (63  -  62) 

Ayr7; 

(c205)  ayr7; 

(d205)  F3  (66  +  65) 

Ayi7; 

(c206)  ayi7; 

(d206) 

F3  (67  -  66) 

Azr7; 

(c207)  azr7; 

(d207)  F3  (69  +  612) 

Azi7; 

(c208)  azi7; 

(d208)  F3  (6)1 -610) 
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5 

I 

I 

After  the  eighth  level  of  computation:  ] 

« 

Axr8;  1 

(c234)  axr8;  < 

l 

I 

I 

(d234)  F3H1  I 

1 

( 

AxiB; 

(c235)  axi8; 

(d235)  F3  H2 

Ayr8; 

(c236)  ayr8; 

(d236)  F3  H3 

Ayl8; 

(c237)  ayi8; 

(d237)  F3  H4 

Azr8; 

(c238)  azr8; 

(d238)  F3  H5 

Azi8; 

(c239J  azi8, 


(d239) 


F3H6 


P.WWPJWJ'P.WWW  WWVWV  V  v  V? 


£ 

i 


0 


V 

V 


After  the  ninth  and  final  level  of  computation: 

Axr9; 

(c265)  axr9; 


(d265)  Jl 


* 

»  * 

A 


!* 


i 


$ 


'v 


Axl9; 

(c266)  axi9; 

(d266)  J2 

Ayr9; 

(c267)  ayr9; 

(d267)  J3 

Ayi9; 

(c268)  ayi9; 

(d268)  •  J4 

Azr9; 

(c269)  azr9; 

(d269)  J5 

Azl9; 

(c27l )  azi9, 

Cd27 1 )  J6 


$ 
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The  text  below  Is  the  actual  commands  used  to  perform  the  calculation  analysis  for  the 
Vector  Potential.  The  commands  are  separated  into  the  computational  levels  developed  in  the 
previous  text. 


These  first  commands  are  in  the  preparation  stage  of  the  process. 

realpart(ax); 

imagpart(ax); 

subst(((((x-xir2M(y-yi)*2)+((z-zi)A2))‘(1/2))lr.JJ); 


Level  one: 

subst(“A1",(x-xi  ),*)$ 
5ubst('A2“.(y-yi),?5)$ 
subst("A3“,(z-zi),£)$ 


Level  two: 

subsU"Bi“."A1“‘2.S)$ 

substf'B2YA2'“2,8)$ 

subst('B3'."A3"2.f?)$ 


Level  three: 

substCC  1  YB  1  ■+'B2"+"B3“,JI); 

Level  Tour: 

subst("Dr,sqrt("CD,R); 

Level  five: 

substCEr,b*"Dr.J?)$ 

substCE2M4*Spi""Dr).f5); 

Level  six: 


substrFi",cosrEr).S)$ 

5ubstrF2'.smrFn,^)$ 

subst('F3M/"E2",!?); 


substCG  1  ",jxri  *  F 1  “,?2)$ 
subst("G2“,jxri*"F2",S)$ 
subst('G3",jxii  *  T  t  ',?*)% 
substCG4",jxii*‘F2\8); 

subst('65“,jyri#  "F  1",S)$ 
substCG6",jyri  *  “F2 ',?»)% 
substCG7",jyii*  "F 1  \S)$ 
substCG8\jyii*‘F2“,8); 

subst(*G9\jzri  * 'F  !".!?)$ 
substCG  10“,Jzri»'F2",^)$ 
substCG  1  r,jzii»’Fr,X)t 
substCG  l2’,Jzii*’F2",X); 


Level  eight: 

substCH  C,"G  I  ■+"G4",J{); 
substCH27'63“-“G2",8); 
substCH3YG5"+'G8'.S); 
substCH4'."G7---G6 
substCH5',“G9“+"G  1 2",!t), 
substCH6',"G1  r-"G10',JJ); 


Level  nine: 

substCJ  I  YF3“*  "H I  ",R); 
substCJ2‘.T3‘«‘H2‘,Jt); 
substCJ3",“F3’*“H3‘,X); 
substCJ4'VF3'*'H4",S), 
substCJ5VF3'*‘H5',8) 
substCJ6YF3"*"H6',X) 
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This  section  illustrates  the  calculation  analysis  of  the  Magnetic  Field. 
The  following  pages  will  show  the  step  by  step  procedure  of  the  computation 
of  the  Parallel  Magnetic  Field. 

The  magnetic  field  in  the  “x"  direction  will  be  the  example  through 
level  one  with  the  complete  complex  six-tuple  shown  for  the  remainder  of 
the  process.  The  text  is  in  MACSYMA  syntax  and  presented  as  the  user  might 
see  it  on  the  screen. 

Hx: 

1  -  %\  b  r 

(-  +  %\  b)  fie  ((jyri  +  Si  jyii)  (z  -  zi)  -  (jzn  +  ?%\  jzii)  (y  -  yi)) 
r 

2 

A  Xpi  r 

rcilpirt(hx)  The  real  part  is  extracted  with  MACSYMA. 

((cos(b  r)  (jyri  (z  -  zi)  -  jzri  (y  -  yi)) 

+  sin(b  r)  (jyii  (z  -  zi)  -  jzii  (y  -  yi)))/r 

-  b  (cos(b  r)  (jyii  (z  -  zi)  -  jzii  (y  -  yi)) 

2 

-  sin(b  r)  (jyri  (z  -  zi)  -  jzri  (y  -  yi))))/(4  Spi  r  ) 

imaflpart(hx)  The  imaginary  part  is  extracted  with  MACSYMA. 

((cos(b  r)  (jyii  (z  -  zi)  -  jzii  (y  -  yi)J 

-  sm(b  r)  (jyri  (z  -  zi)  -  jzri  (y  -  yi)))/r 
b  tcos(b  r)  (jyri  (z  -  zi)  -  jzri  (y  -  yi)) 

o 

im 

*  sin(b  r)  (jyii  (z  -  zi)  -  jzii  (y  -  yi))))/(4  Rpi  r  ) 
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Here  we  begin  the  computation  process 

The  example  below  recalls  the  real  part  of  Hx,  then  substitutes  the  expanded  version  of  "r"  into  the 
expression. 

rcilpirt(hx); 

hxr; 

((cos(b  r)  (jyri  (z  -  zi)  -  jzri  (y  -  yiJ) 

+  sin(b  r)  (jyii  (z  -  zl)  -  jzii  (y  -  yi)))/r 
-  b  (cos(b  r)  (jyii  (z  -  zi)  -  jzii  (y  -  yi )) 


n 

A. 

-  sinib  r)  (jyri  (z  -  zi)  -  jzri  (y  -  yi))))/(4  Xpi  r  ) 


3ubst(((((x->dr2)+((y-ylF2)+((z-zi)*2))*(1/2)),r,%); 

HxrO; 


t  ■ 


> 


2  2  2 

((cos(b  sqrti(z  -  zi)  +  (y  -  yi)  +  (x  -  xi)  )) 

(jyri  (z  -  zi)  -  jzri  (y  -  yi))  +  sin(b 
2  2  2 

sqrt((z  -  Zl)  +  (y  -  yi)  +  (x  -  xi)  ))  (jyn  (z  -  zi)  -  jzii  (y  -  yi))) 

2  2  2 
/sqrt((z  -  zi)  +  (y  -  yi)  +  (x  -  xi)  ) 

2  2  2 

-  b  (cos(b  sqrt((z  -  zi)  +  (y  -  yi)  ♦  (x  -  xi)  )) 

(jyii  (z  -  zi)  -  jzii  (y  -  yi))  -  sm(b 
2  2  2 

sqrt((z  -  zi)  *  (y  -  yi)  +  (x  -  xi)  ))  (Jyri  (z  -  zi)  -  jzri  (y  -  yi)))) 

2  2  2 
/(4  3pi  ((z  -  zi)  ♦  (y  -  yi)  +  (x  -  xi)  )) 


% 

% 


V, 

V 
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After  the  first  level  of  computation  the  real  part  of  the  X 
directed  Magnetic  Field  has  the  form 


hxrl; 


2  2  2 


(C(jyii  A3  -  jzii  A2)  sin(b  sqrt(A3  +  A2  +  A1  )) 

2  2  2  2  2  2 
+  (jyri  A3  -  jzri  A2)  cos(b  sqrt(A3  +  A2  +  A1  )))/sqrt(A3  +  A2  +  A1  ) 

2  2  2 

-  b  ((jyil  A3  -  jzii  A2)  cos(b  sqrt(A3  +  A2  +  A1  )) 

2  2  2  2  2  2 

-  (jyri  A3  -  jzri  A2)sin(b  sqrt(A3  +  A2  +  A1  ))))/(4  Spi  (A3  +  A2  +  A I  )) 


After  the  second  level  of  computation: 

hxr2; 

IB 1 5  -  B20 )  sin(b  sqrt(B3  +  B2  +  BD)  +  (812  -  B17)  cos(b  sqrt(B3  ♦  B2  ♦  B 1 )) 

( - 

sqrt(B3  +  B2  >  B 1 ) 

-  b  ((B 1 5  -  B20)  cos(b  sqrt(B3  >  B2  ♦  BD) 

-  (B 1 2  -  B 1 7)  sin(b  sqrt(B3  *  B2  ♦  B 1  ))))/{4  Rpi  (B3  +  B2  +  B 1 )) 

hx  1 2; 

(b  (CB 1 5  -  B20)  sin(b  sqrt(B3  ♦  B2  +  B D) 

+  (B12  -  B17)  cos(b  sqrt(B3  ♦  B2  +  BD)) 

+  «BI5  -  B20)  cos(b  sqrt(B3  ♦  B2  ♦  B D) 

-  (B 12  -  B 1 7)  sm(b  sqr((B3  +  62  +  B I  ))Vsqrt(B3  +  B2  +  B  D) 


/(4  Spi  (B3  ♦  B2  ♦  BD) 


hyr2; 

sm(b  sqrt(B3  ♦  B2  +  BUHB19  -  B9)  +  cos(b  sqrt(B3  +  B2  +  BD)  (B16  -  B6) 


sqrt(B3  ♦  62  «•  B 1 ) 

-  b  (cos(b  sqrt(B3  ♦  B2  ♦  BD)  (B19  -  B9) 

-  5in(b  sqrt(B3  ♦  B2  ♦  B 1 ))  (B 1 6  -  B6  )))/(4  Xpi  (B3  ♦  B2  ♦  B I )) 


hy12; 

(b  (sin(b  sqrt(B3  ♦  B2  ♦  B 1 ))  (B 1 9  -  B9) 

♦  cos(b  $qrt(B3  +  B2  +  B 1 ))  (B 16  -  B6)) 

+  (cos(b  sqrt(B3  +  B2  +  B 1 ))  (BI9  -  B9) 

-  sin(b  sqrt(B3  ♦  B2  ♦  BD)  (B16  -  B6))/sqrl(B3  +  B2  +  B  D) 
/{A  Spi  (B3  +  B2  +  BD) 


hzr2; 

sm(b  sqrt(B3  ♦  B2  ♦  B 1 ))  (B8  -  B 1 3)  ♦  cos(b  sqrt(B3  ♦  B2  ♦  B 1 ))  (B5  -  B 1 0) 
sqrt(B3  +  B2  +  B 1 ) 

-  b  (cos(b  sqrt(B3  ♦  B2  ♦  81 ))  (B8  -  B13) 

-  sin(b  sqrt(B3  +  B2  +  BD)  (85  -  B10JM4  Spi  (B3  ♦  B2  ♦  BD) 


hzi2 

(b  (sin(b  sqrt(B3  +  B2  B 1 ))  (B8  -  B 1 3) 

+  cos(b  sqrt(B3  +  B2  +  BD)  (B5  -  BIO)) 

♦  (cos(b  sqrt(B3  ♦  B2  ♦  B1 ))  (B8  -  B13) 

-  sin(b  sqrt(B3  ♦  B2  +  B 1))  (B5  -  B 1 0))/sqrt(B3  +  B2  +  B I )) 


After  the  third  level  of  computation: 

hxr3; 

cos(b  sqrKC  1 ))  C3  +  sin(b  sqrt(C  1 ))  C2 


sqrt(C  1 ) 

-  b  (cos(b  sqrt(C  1 ))  C2  -  sin(b  sqrt(C  I ))  C3))/(4  Rpi  C 1 ) 

hxi3; 

cos(b  sqrt(C  1 ))  C2  -  sin(b  sqrt(C  1 ))  C3 
sqrKC  1 ) 

+  b  (cos(b  sqrt(C  1 ))  C3  *  sln(b  sqrKC ! ))  C2))/(4  3pi  C 1 ) 


hyr3; 

coslb  sqrt(C  I ))  C5  +  sin(b  sqrt(C  1 ))  C4 

( - 

sqrKC  1 ) 

-  b  (cos(b  sqrt(C  1 ))  C4  -  sin(b  sqrt(C  1 ))  C5))/(4  Spi  C 1 ) 

hyi3; 

cos(b  sqrt(C  1 ))  C4  -  sin(b  sqrt(C  1 ))  C5 

( - 

sqrt(C  I ) 

+  b  (cos(b  sqrt(C  1 ))  C5  +  sin(b  sqrKC  1 ))  C4M4  J5pl  C 1 ) 


hzr3; 

cos(b  sqrt(C  1 ))  C7  ♦  sin(b  sqrt(C ! ))  C6 
sqrt(C  i ) 

-  b  (cos(b  sqrt(C  1 ))  C6  -  sin(b  sqrKC  1 ))  C7M4  ?5pi  C 1 ) 

hzl3; 

cos(b  sqrt(C  1 ))  C6  -  sm(b  sqrt(C  I ))  C7 

( - 

sqrt(C  1 ) 

+  b  (cos(b  sqrt(C  I ))  C7  +  sinCb  sqrtCC  1 ))  C6))/(4  ?5pj  C 1 ) 
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D2 


hzr4; 


C6  sin(b  Dl)  +  C7  cos(b  D!) 


After  the  fifth  level  of  computation: 

hxr5; 

((C2  sin(E  I )  +  C3  cos(E  1 ))  E2  -  b  (C2  cos(E  1 )  -  C3  stn(E  1 )))  E3 

hxiS; 

((C2  cos(E  1 )  -  C3  sin(E  1 ))  E2  +  b  (C2  sin(E  1 )  +  C3  cos(E  1 )))  E3 

hyr5: 

((C4  sin(E I )  ♦  C5  cos(E I ))  E2  -  b  (C4  cos(£ I)  -  C5  sin(E  1 )))  E3 

hyl5; 

((C4  cos(E  I )  -  C5  sin(E  D)  E2  +  b  (C4  sin(E  1 )  +  C5  cos(E  1 )))  E3 

kzrS; 

((C6  sin(E I)  +  C7  cos(E  1 ))  E2  -  b  (C6  cos(E  1 )  -  C7  3in(E I )))  E3 

hzl5; 

((C6  cos(E  1 )  -  C7  sin(E  1 ))  E2  ♦  b  (C6  sin(E  1 )  ♦  C7  cos(E  1 )))  E3 

After  the  sixth  level  of  computation: 

hxr6; 

E3  (E2  (C2  F2  +  C3  FI)  -  b  (C2  FI  -C3F2)) 

hxi6; 

E3(E2  (C2FI  -  C3  F2)+b(C2  F2  +  C3FD) 

hyr6; 

E3  (E2(C4  F2  +  C5  FI)  -  b  (C4  FI  -C5F2)) 

hyi6; 

E3(E2  (C4FI  -  C5  F2)+b(C4  F2  +  C5FU) 

hzr6; 

E3(E2  (C6  F2  +  C7  F ! )  -  b  (C6  F 1  -C7F2)) 

hzi6; 


E3(E2  (C6FI  -  C7  F2)  +  b  (C6  F2  +  C7FD) 


After  the  seventh  level  of  computation: 

hxr7; 

E3  (E2  (67  *  62)  -  b  (61  -68)) 

hxi7; 

E3  (E2C61  -  68)  ♦  b  (67  +  62)) 

hyr7; 

E3  (E2  (69  ♦  64)  -b  (63  -610)) 

hyi7; 

E3  (b  (69  >  64)  t  E2  (63  -  610)) 

hzr7; 

E3  (E2  (66  +  61 1 )  -  b  (65  -  612)) 

hzl7; 

E3  (b  (66  +  61D  +  E2  (65-612)) 

After  the  eighth  level  of  computation: 

hxr8; 

E3(E2  HI  -  bH4) 

hxi8; 

E3  (E2  H4  +  b  H 1 ) 

hyr8; 

E3  (E2  H2  -  b  H5) 

hyi8; 

E3  (£2  H5  +  b  H2) 

hzr8; 

E3  (E2  H3  -  b  H6) 

hzl8; 

E3  (E2  H6  ♦  b  H3) 
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After  the  ninth  level  of  computation: 

hxr9; 


hx!9; 


hyr9; 


hyi9; 


hzr9; 


E3  (J1  -J10) 


E3  (J7  +  J4) 


E3  (J2-J1I; 


E3  (J8  f  J5) 


E3(J3-JI2) 


hz!9; 

E3  (J9  -►  J6) 

After  the  tenth  level  of  computation: 


hxrlO; 

hxilO; 

hyrtO; 

hyilO; 

hzriO; 

hzilO: 


E3KI 

E3K4 

E3K2 

E3K5 

E3K3 

E3K6 
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fcU.  *<■*  lU'lC 4l.  it  St«. 


The  text  below  is  ihe  actual  commands  used  to  perform  the  calculation  analysis  for  the 
Magnetic  Field.  The  commands  are  separated  into  the  computational  levels  developed  in  the 
previous  text. 

These  first  commands  are  in  the  preparation  stage  of  ihe  process. 

realparUhx); 

imagpart(hx); 

subst(((((x-xi  )‘2  W(y-yi  )‘2)+((z-zt  T2  ))*(  1  /2)).r  ,S); 


Level  one: 

subst('Ar.(x-xi),X); 

substCA2",(y-yi),J5); 

substCA3",(z-z1),X); 


Level  two 

i 

substCB  1", "A  r*2.S); 
subst('B2’.'A2'*2.JJ); 
subst(“B3"."A3“2.R); 

SubstCB4",jxri*  "A  1  ".X); 

not  used 

V\ 

substCB5',jxri*  “A2".8); 

a 

•> 

subst("B6\jxri*  "A3".K); 

^  ■ 

am 

subst('B7'.jxii"‘Ar,R), 
su0stC"B8'.jxii»  "A2".?S); 

not  used 

? 

subsU"B9",jxii*  "A3',R); 
substCB  lO'.jyri*  'A  l“,%); 
substCB  1  f.jyri*  "A2",/S), 

not  used 

/* 

subst("B12".jyri*'A3',X); 
subst("B  1 3",jyii*  "A  I'.S), 
substCB  1 4',jyii«  'A2‘,X); 
substCB  15“,jyii*  “A3".K); 

not  used 

►i 

3! 

substCB  1 6",  jzri  *  "A  1  ",X); 
substCB  1 7",jzri*  ’A2“.?J), 
substCB  1 8". jzri  *  "A3", S); 
subst("B  1 9",jzti*  "A  1  ",R); 

not  used 

substCB20‘.jzii*  _A2'.S); 
substCB2r,jzii"‘A3'.J5). 

not  used 

■2 

i 


This  section  illustrates  the  calculation  analysis  of  the  Electric  Field. 

The  following  pages  will  show  the  step  hy  step  procedure  of  the  computation 
of  the  Parallel  Electric  Field. 

The  Electric  Field  in  the  "x"  direction  will  be  the  example  through  level 
four  with  the  complete  complex  six-tuple  shown  for  the  remainder  of  the 
process.  The  text  is  in  MACSYMA  syntax  and  presented  as  the  user  might  see 
it  on  the  screen. 


Ex; 

(c9)  ex; 

-Sibr  3  %i  b  1  2  2  2 

(d9)  Se  ((jxri  +  Si  jxii)  ( - +  b  )  ((2  -  zi)  +  (y  -  yi) ) 

r  2 
r 

3  Si  b  3  2 

+  ( - +  —  b  )  ((jzri  +  Si  jzii)  (x  -  xi)  (z  -  zi) 

r  2 
r 


3 

+  (jyri  +  Si  jyii)  (x  -  xi)  (y  -  yi)))/(4  Spi  r  ) 


realaaruex)  The  real  part  is  extracted  with  MACSYtlA. 
(clO)  exr; 

2  1  3  b  jxri  2  2 

(d  10)  (sin(h  r)  ((Jxii  (b  — ) - )  ((z  -  zi)  +  (y  -  yi) ) 

2  r 


3  b  (jzrl  (x  -  xi)  (z  -  zl)  +  Jyrl  (x  -  xl)  (y  -  yi)) 


3  2 

+  ( —  b  )  (jzll  (x  -  xl)  (z  -  zl)  +  Jyll  (x  -  xl)  (y  -  yi))) 
2 


3  b  jxii  2  1  2  2 

+  cos(b  r)  (( - ♦  jxri  (b  — ))  ((z  -  zi)  ♦  (y  -  yi) ) 

r  2 


3  2 

+  ( - b  )  (jzri  (x  -  xi)  (z  -  zi)  +  jyri  (x  -  xi)  (y  -  yi)) 

2 


3  b  (jzii  (x  -  xi)  (z  -  zi)  +  jyii  (x  -  xi)  (y  -  yi)) 


-))/(4  3pi  r  ) 


imagpart(ex)  The  imaginary  part  is  extracted  with  MACSYMA. 
(ell)  exi; 

2  1  3  b  jxri  2  2 

(d!  1)  (cos(b  r)  ((jxii  (b - ) - )((z-zt)  +  (y  -  yi) ) 

2  r 


3  b  (jzri  (x  -  xi)  (z  -  zi)  +  jyri  (x  -  xi)  (y  -  yi)) 


♦  ( - b  )  (jzii  (x  -  xi)  (z  -  zl)  +  jyii  (x  -  xl)  (y  -  yi))) 


3  b  jxii  2  1  2  2 

-  sin(b  r)  (( - ♦  jxri  (b  — ))  ((z  -  zi)  +  (y  -  yi)  ) 

r  2 


3  2 

♦  ( —  b  )  (jzri  (x  -  xi)  (z  -  zi)  +  jyri  (x  -  xi)  (y  -  yi)) 
2 


3  b  (jzii  (x  -  xi)  (z  -  zi)  ♦  jyii  (x  -  xi)  (y  -  yi)) 


-M4  Xpi  r  ) 


Here  we  begin  the  comPuUUQn  process: 

The  example  below  recalls  the  real  part  of  Ex,  then  substitutes  the  expanded  version  of  "r"  into  the 
expression. 


I 

f  . 


ExrO;  Substitute  the  expanded  version  of  "r‘  into  the  expression, 
(c  1 2)  exrO; 

2  2  2 

(d  12)  (sm(b  sqrt((z  -  zi)  +  (y  -  yi)  +  (x  -  xi)  )) 

3  b  (jzri  (x  -  xi)  (z  -  zi)  +  jyri  (x  -  xi)  (y  -  yi)) 

2  2  2 
sgrt((z  -  zi)  ♦  (y  -  yi)  ♦  (x  -  xi)  ) 

2  1 

J  2  2  2 
(z  -  zi)  +  (y  -  yi)  ♦  (x  -  xi) 

3  b  jxri  2  2 

- )((z-zi)  +  (y  -  yi)  ) 

2  2  2 
sqrt((z  -  zi)  +  (y  -  yi)  +  (x  -  xi) ) 

3  2 

+  ( - b  ) 

2  2  2 
(z  -  zi)  +  (y  -  yi)  +  (x  -  xi) 


(jzii  (x  -  xi)  (z  -  zi)  -*•  jyli  (x  -  xi)  (y  -  yi))) 

2  2  2 

+  cos(b  sqrt((z  -  zi)  ♦  (y  -  yi)  ♦  (k  -  xi)  )) 

3  b  (jzii  (x  -  xi)  (z  -  zi)  ♦  jyii  (x  -  xi)  (y  -  yi)) 

2  2  2 
sqrt((z  -  zi)  +  (y  -  yi)  +  (x  -  xi)  ) 

3  b  jxii 

♦  ( - 

2  2  2 
sqrt((z  -  zi)  ♦  (y  -  yi)  +  (x  -  xi)  ) 


2  1  2  2 

♦  jxri  (b  - ))  ((z  -  zi)  +  (y  -  yi)  ) 

2  2  2 
(z  -  zi)  +  (y  -  yi)  ♦  (x  -  xi) 

3  2 

♦  ( - b  ) 

2  2  2 


(z  -  zi)  ♦  (y  -  yi)  +  (x  -  xi) 

2  2  23/2 

(jzri  (x  -  xi)  (z  -  zi)  +  jyri  (x  -  xi)  (y  -  yi))))  /(4  Spi  ((z  -  zi)  ♦  (y  -  yi)  ♦  (x  -  xi)  )  ) 


V 

v 
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ExIO;  Substitute  the  expanded  version  of  "r“  into  the  expression, 
(c  13)  exiO, 

2  2  2 

(d  1 3)  (cos(b  sqrt((z  -  zi)  +  (y  -  yi)  +  (x  -  xi) )) 

3  b  (jzri  (x  -  xi)  (z  -  zi)  +  jyri  (x  -  xi)  (y  -  yi)) 

( - 

2  2  2 
sqrt((z  -  zi)  +  (y  -  yi)  ♦  (x  -  xi) ) 

2  1 

]  2  2  2 
(z  -  zi)  +  (y  -  yi)  +  (x  -  xi) 

3  b  jxri  2  2 

- )  ((2  -  2j)  +  (y  -  yj)  ) 

2  2  2 
sqrt((z  -  zi)  +  (y  -  yi)  +  (x  -  xi)  ) 

3  2 

+  ( - b  ) 

2  2  2 
(z  -  zi)  +  (y  -  yi)  +  (x  -  xi) 

(jzii  (x  -  xi)  (z  -  zi)  +  jyii  (x  -  xi)  (y  -  yi))) 

2  2  2 

-  sin(b  sqrtttz  -  zi)  +  (y  -  yi)  +  (x  -  xi)  )) 

3  b  (jzii  (x  -  xi)  (z  -  zi)  ♦  jyii  (x  -  xi)  (y  -  yi)) 

2  2  2 
sqrt((z  -  zi)  +  (y  -  yi)  +  (x  -  xi)  ) 

3  b  jxii 


2  2  2 
sqrt((z  -  zi)  +  (y  -  yi)  ♦  (x  -  xi) ) 


2  1  2  2 

*  jxri  (b  - ))((z-zi)  +(y-yi)  ) 

2  2  2 
(z  -  zt)  +  (y  -  yi)  ♦  (x  -  xi) 

3  2 

♦  ( - b  ) 

2  2  2 
(z  -  zi)  +  (y  -  yi)  +  (x  -  xi) 

(Jzri  (x  -  xi)  (z  -  zi)  +  jyri  (x  -  xi)  (y  -  yi)))) 

2  2  2  3/2 

/(4  JJpi  ((z  -  zi)  +  (y  -  yi)  *  (x  -  xi)  )  ) 


After  the  first  level  of  computation: 


Exrl; 

lc14)  exrl; 


2  2  2  1  3  b  jxri 

Id  141  (((A3  +  A2  )  (jxii  (b - ) - ) 

2  2  2  2  2  2 


A3  +  A2  +  A1  sqrt(A3  +  A2  +  A1  ) 
3  2 

+  (jzii  A 1  A3  +  jyii  A 1  A2)  ( - - b  ) 

2  2  2 
A3  +A2  ♦  A1 


3b(jzri  A1  A3  +  jyri  A I  A2)  2  2  2 

+ - )  sin(b  Sqrt(A3  +  A2  ♦  A 1  )) 

2  2  2 
sqrt(A3  +  A2  ♦  A 1  ) 


2  2  2  1  3  b  jxii 

♦  ((A3  +  A2  )  (jxri  (b -  )  + 

2  2  2  2  2  2 
A3  ♦  A2  +  A 1  sqrt(A3  +  A2  +  A 1  ) 

3  2 

♦  (jzri  A1  A3  +  jyri  A1  A2)  ( - b  ) 

2  2  2 
A3  ♦  A2  +  A 1 

3  b  (jzii  A 1  A3  +  jyii  A 1  A2)  2  2  2 

- )  cos(t)  5qrt(A3  +  A2  +  A I  ))) 

2  2  2 
sqrt(A3  +  A2  +  A 1  ) 

2  223/2 

A 4  Spi  (A3  +  A2  ♦  A1  )  ) 
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(cl5)  exi  1 ; 


2  2  2  1  3  b  jxri 

Cd  15)  (((A3  +  A2  )  (jxii  (b - ) - 

2  2  2  2  2  2 
A3  +  A2  +  A 1  sqrt(A3  +  A2  +  A 1  ) 

3  2 

+  (jzii  A1  A3  +  jyii  A1  A2)  ( - b  ) 

2  2  2 
A3  +  A2  +  A! 

3  b  (jzri  A1  A3  ♦  jyri  A1  A2)  2  2  2 

+ - )  C0S(b  sqrt(A3  +  A2  +  A 1  )) 

2  2  2 
sqrt(A3  +  A2  +  A I  ) 

2  2  2  I  3  b  jxii 

-((A3  +  A2  )(jxri(b - )■*■ - ) 

2  2  2  2  2  2 
A3  +  A2  ♦  A1  sqrt(A3  +A2  +  A1  ) 

3  2 

♦  (jzri  A 1  A3  +  jyri  A 1  A2)  ( - b  ) 

2  2  2 
A3  +A2  +  A1 

3  b  (jzii  A1  A3 +'jyii  Ai  A2)  2  2  2 

- )  sin(b  sqrt(A3  +  A2  +A1  ))) 

2  2  2 
sqrt(A3  >A2  -*-A1  ) 


2  223/2 

/(4  Kpi  (A3  ♦  A2  +  A)  )  ) 


After  the  second  level  of  computation: 

Exr2; 

(cl 8)  exr2; 

2  1  B22 

(d  18)  (((83  +  82)  (jxii  (b - ) - ) 

B3  +  B2  +  B1  sqrt(B3  +  B2  ♦  B 1 ) 

3  2  3  b  (A3  B 16  +  A2BJ0) 

+  (A3B19  + A2B13K - b)  + - ) 

B3  +  B2  +  BI  sqrt(B3  +  B2  +  B 1 ) 

sin(b  sqrt(B3  +  B2  ♦  B 1 ))  +  ((B3  +  B2) 

2  1  B23 

(jxri  (b - )  + - ) 

B3  +  B2  +  B1  sqrt(B3  +  B2  +  B ) ) 

3  2  3b(A3Bl9  + A2B13) 

+  (A3B16  ♦  A2B10M - b  ) - ) 

B3  +  B2  +  B!  sqrt(83  +  B2  +  B 1 ) 

3/2 

cos(b  sqrt(B3  +  B2  ♦  B 1  )))/(4  8pi  (B3  ♦  82  ♦  B 1 )  ) 


Exl2; 

(c  19)  exi2; 


2  1  B22 

(d19)(((B3  +  B2)  (jxii  (b - ) - ) 

B3  +  B2  +  B!  sqrt(B3  ♦  B2  +  B 1 ) 

3  2  3  b  (A3  B16  +  A2  BIO) 

♦  (A3B19  +  A2B13M - b  )♦ - ) 

B3  +  B2  +  B1  sqrt(B3  +  B2  +  B 1 ) 

cos(b  sqrt(B3  +  B2  +  B 1 ))  -  ((63  +  B2) 

2  1  B23 

(jxri  (b - )  + - ) 

B3  +  B2  +  B1  sqrt(B3  +  B2  +  B 1 ) 

3  2  3  b  (A3  BI9  +  A2  B13) 

♦  (A3BI6  +  A2BI0M - b  ) - ) 

B3  +  B2  +  B1  sqrt(B3  ♦  B2  B1) 


3/2 

sln(b  sqrt(83  ♦  B2  ♦  B 1  )))/(4  Rp(  CB3  •*■62  +  61)  ) 


a 


R 


After  the  third  level  of  computation: 


Exr3; 

(c22)  exr3; 


* 


s  ' 


•r. 

A 


s 


s 


■i 


-3 


3  b (C38  +  C36)  3  2 

(d22)  (sln(b  sqrt(C  1 ))  ( - +  (-  -  b  )  (C37  +  C35) 

sqrt(CI)  Cl 

2  1  B22 

♦  (jxii  (b  —  ) - )  C33)  +  costb  sqrt(C  1 )) 

Cl  sqrt(CI) 

3  2  3  b  (C37  *  C35)  B23  2  1 

((--  -  b  )  (C38  ♦  C36) - +  ( - +  jxri  (b  -  -  ))  C33)) 

Cl  sqrt(CI)  sqrt(Cl)  Cl 

3/2 

/(4SpiC1  ) 


Exl3; 

(c23)  exi3, 


3  b  (C38  ♦  C36)  3  2 

(d23)  (cos(b  sqrt(C  1 ))  ( - +  (  --  -  b  )  (C37  ♦  C35) 

sqrt(C1)  Cl 

2  I  B22 

+  (jxii  (b  —  ) - )  C33)  -  sin(b  sqrt(C  1 )) 

C 1  sqrt(C  1 ) 

3  2  3  b  (C37  +  C35)  B23  2  ! 

((--  -  b  )  (C38  ♦  C36) - - ♦  ( - ♦  jxri  (b  -  -  ))  C33)) 

Cl  sqrt(CI)  sqrKC  \ )  Cl 

3/2 

/(4F.pi  Cl  ) 
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After  the  forth  level  of  computation: 


Exr4; 

(c51)  exr4; 

2  3  b  D34  2  B22 

(d5 1 )  (sin(b  0 1 )  ((3  D33  -  b  )  035  ♦ - +  C33  (jxii  (b  -  D33) - )) 

D1  D1 


3  b  D35  2  2  B23 

+  cos(b  D 1 )  ( - +  (3  033  -  b  )  D34  +  C33  (jxri  (b  -  033)  +  —  ))) 

D1  D1 


✓(01  D2) 

Exl4; 

(c52)  exi4; 


2  3  b  D34  2  B22 

(d52)  (cos(b  D 1 )  ((3  D33  -  b  )  035  + - +  C33  (jxii  (b  -  D33) - )) 

D1  D1 

3  b  035  2  2  B23 

-  sin(b  DIM - +  (3  D33  -  b  ) D34  +  C33  (jxri  <b  -D33)  + —  ))) 

01  01 

/(D1  02) 

After  the  fifth  level  of  computation: 


Exr5; 

(c53)  exr5, 


2 

(dS3)  E2  E3  (cos(E  1 )  (-  E2  E7  +  034  (E5  -  b  )  +  C33  (jxri  E4  +  B23  E2)) 

2 

♦  s!n(E  1 )  (E2  E6  +  035  (E5  -  b  )  +  C33  (jxii  E4  -  B22  E2)» 

Exl5; 

(c54)  exi5, 

2 

(d54)  E2  E3  (cos(E  1 )  (E2  E6  ♦  035  (E5  -  b  )  ♦  C33  (jxii  E4  -  B22  E2)) 

2 

-  sm(E  1)  (-  E2  E7  +  034  (E5  -  b  )  ♦  C33  (jxri  E4  ♦  B23  E2)» 
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Eyr5; 

(c55)  eyr5; 


2 

(055)  E2  E3  (cos(E  1 )  (-  E2  E9  +  D36  (E5  -  b  )  +  C32  (jyrl  E 4  +  B25  E2)) 

2 

+  sin(E  1 )  (E2  E8  +  037  (E5  -  b  )  ♦  C32  (jyii  E4  -  B24  E2))) 


EyJ5; 

(c57)  eyi5; 


2 

(057)  E2  E3  (cos(E  I )  (E2  E8  ♦  037  (E5  -  5  )  ♦  C32  (jyii  E4  -  B24  E2)) 

2 

-  sin(E  1 )  (-  E2  E9  +  036  (E5  -  b  )  +  C32  (jyri  E4  +  B25  E2))) 


EzrS; 

(c58)  ezr5, 


2 

(058)  E2  E3  (sm(E  1 )  (039  (E5  -  b  )♦  C3 1  (jzii  E4  -  B26  E2)  +  E10  E2) 

2 

*  cos(E  1 )  (D38  (E5  -  b  )  ♦  C3 1  (jzri  E4  ♦  B27  E2)  -Ell  E2)) 


Ezl5; 

(c59)  ezi5. 


(059 )  E2  E3  (costE  I )  (D39  ^E5  -  b  )  ♦  C31  (jzii  E4  -  B26  E2)  ♦  E 10  E2) 

i. 

-  sm(E  l )  (D36  (E5  -  b  )*  C3 »  (jzri  E4  ♦  B27  E2)  -  E II  E2)) 


After  the  sixth  level  of  computation: 

Exr6; 

(c60)  exr6; 

(d60)  F3  (F2  (C33  (F0  -  F5)  +  035  FA  +  E6)  +  FI  (C33  (F7  +  F6)  ♦  D34  FA  -  E7)) 

Exl6; 

(c61)  exi6; 

(d6 1 )  F3  (F I  (C33  (F8  -  F5)  +  035  F4  ♦  E6)  -  F2  (C33  (F7  +  F6)  +  D34  F4  -  E7)) 

Eyr6; 

(c62)  eyr6; 

Cd62)  F3  (F2  (C32  (F 12  -  F9)  +  D37  F4  ♦  E8)  ♦  F 1  (036  F4  ♦  C32  (F 11  ♦  F 10)  -  E9)) 

Eyl6; 

(c63)  eyi6; 

(d63)  F3  (F I  (C32  (F 12  -  F9)  +  037  F4  +  E8)  -  F2  (D36  F4  +  C32  (F 1 1  ♦  F 10)  -  E9)) 

Ezr6; 

(c64)  ezr 6; 

(d64)F3  (F2  (D39F4  +  C31  (F 16  -  F 13)  +  E 10)  +  F 1  (D38F4  +  C31  (F15  ♦  F 14)  -  E 1 1 )) 

Ezl6; 

(c65)  e’i6; 

ld65)  F3  (F I  (039F4  +  C3I  (F 16  -  F 13)  +  E 10)  -  F2  (D38  F4  +  C31  (F 15  ♦  F 14)  -  E 11 » 


After  the  seventh  level  of  computation: 


Exr7; 

(c66)  exr7; 

(<166 )  F3  (F 1  (622  ♦  C33  G2 1  -  E7)  +  F2  (623  ♦  C33  620  ♦  E6H 

Exl7; 

(c67)  exi7; 

(d67)  F3  (F 1  (623  ♦  C33  620  ♦  E6)  -  F2  (622  +  C33  62 1  -  E7M 

Eyr7; 

(c68)  eyr7; 

(468)  F3  (F 1  (632  +  C32  63 1  -  E9)  ♦  F2  (633  +  C32  G30  ♦  E8)) 

Eyl7; 

(c69)  eyi7| 

(469)  F3  (F 1  (633  ♦  C32  G30  +  E8)  -  F2  (632  ♦  C32  63 1  -  E9» 

Ezr7; 

(c70)  ezr7; 

(470)  F3  (FI  (G42  +  C3I  641  -  Ell)*  F2  (643  ♦  C3 1  640  ♦  E 10)) 

Ezl7; 

(c7l )  czi7; 

(471)  F3  (FI  (643  +  C31  640  +  EI0)-F2(642  ♦  C31  641  -  Ell)) 


After  the  eighth  level  of  computation: 

Exr8; 

(c72)  exr8; 

(472)  F3  (F2  (H23  +  H21 )  ♦  FI  (H22  +  H20)) 

Exi8; 

(c73)  exi8; 

(473)  F3  (F 1  (H23  +  H21)-F2(H22  ♦  H20)) 

Eyr8; 

(c 74)  eyr8, 


(474) 


F3  (F2  (H33  +  H3 D  +  F 1  (H32  +  H30)) 


WWWWWifWWiWW 


Ift 

L 


•  w 

.-> 

Eyl8; 

(c75)  eyi8; 

1! 

(d75) 

F3  (F 1  (H33  ♦  H3 1 )  -  F2  (H32  +  H30)) 

Ezr8; 

£ 

(c76)  ezr8. 

V 

(d76) 

F3  (F2  (H43  ♦  H41 ) ♦ FI  (H42  ♦  H40)) 

« 

.•* 

Ezl8; 

(c77)  e?i8. 

i? 

<<J77> 

C3  (f  1  (H43  ♦  H41 )  -  F2  (H42  ♦  H40)) 

V 

After  the 

ninth  level  of  computation: 

V 

■/ 

E*r9; 

(c78)  exr9. 

(d78) 

F3 (F 1  J21  +  F2  J20) 

Exi9; 

(c79)  exi9. 

■ 

(d79) 

F3  (FI  J20-F2J2I) 

Eyr9; 

V 

(c80 1  »yr9. 

(d80) 

F3  (F 1  J31  ♦  F2J30) 

■ 

Eyi9; 

(cQ  1 1  eyi9; 

S' 

> 

ft 

(d8U 

F3  (F 1  J30-F2J3I) 

Ezr9; 

v» 

(c82)  e.’p9. 

\' 

(d82) 

F 3  (F  I  J41  +  F2  J40) 

s 

»' 

Ez«9; 

(c83>  ezi9. 

(d83) 

F3(F!  J40-F2J41) 

'  ft 


> 

V, 
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After  the  tenth  level  of  computation: 

ExrlO; 

(c84)  exrlO; 

(d84)  F3  (K21  ♦  K20) 

Exi 10; 

(c85)  exi  10; 

(d85)  F3  (K22-K23) 

EyrIO; 

(c86)  eyrIO; 

(d86)  F3  (K31  ♦  K30) 

EyilO; 

(c87)  eyi  10. 

(d87)  F3  (K32  -  K33) 

EzrlO: 

(c88)  ezrlO; 

(d88)  F3  (KC4 1  +  K40) 

EzitO; 

(c89)  ezi  10; 

(d89)  F3  (K42  -  K43) 

After  the  eleventh  level  of  computation: 

Exr  1 1 ; 

(c90)  exrl  1; 

(d90)  F3  L20 

Ext  f 1; 

(c91 )  exit  1; 

ld91)  F3L21 

Eyrl 1; 

(c92)  eyrl  1, 

(d92)  F3L30 
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iU-t-.4'«  »«*  *  *^»  *A  X.t  i A  *.*  A»V  «■*. 


Eyl  1 1 : 

(c93)  eyill; 


Ezr 1 1 ; 

(c94)  ezr  1 1 ; 


Ezi  1 1 ; 

(c95)  ezi  11, 


F3L31 


F3L40 


(d95)  F3L4t 

After  the  twelfth  and  final  level  of  computation: 

Exr12; 

(c96)  exr12; 


Exl  12: 

(c97)  exi  12, 


Eyr 12; 

(c98)  eyr  12; 


EyM2; 

(c99)  eyi  12; 


Ezr 12; 

(c  1 00)  ezr  12; 

(d  100) 

Ezi  1 2: 

(clOD  ezi  12. 

Cd  10 1 ) 


Ex  real 


Ex  imag 


Ey  real 


Ey  Imag 


£z  real 


Ez  imag 


^  V 


!W5 


MACSYMA  Command  Lifting  for  EgmUfij  Electric  Field 

The  text  below  is  the  actual  commands  used  to  perform  the  calculation  analysis  for  the 
Electric  Field  The  commands  are  separated  into  the  computational  levels  developed  in  the 
previous  text. 

Note  here  that  certain  MACSYMA  conventions  require  what  would  seem  a  very  round 
about  way  to  substitute  into  equcfons.  For  example,  in  level  two.  notice  the  required  steps  to 
perform  the  operations  that  were  quite  straight  forward  in  the  Magnetic  field  computations. 
These  extra  steps  were  necessary  because  of  the  way  these  equations  were  loaded  into 
MACSYMA  originally.  The  only  impact  on  ihe  behavior  of  the  calculation  flow  is  that  certain 
choices  of  variable  combinations  are  made  by  the  order  of  these  input  commands  and  as  such, 
other  combinations  are  certainly  possible  and  indeed  may  have  been  more  natual.  However, 
the  overall  behavior  of  the  equations  proved  to  be  symmetric  and  therefore  caused  no  real 
concern. 

The  substitution  commands  that  were  identical  to  those  of  ihe  Magnetic  Field 
commands  are  in  standard  print.  The  commands  required  that  were  strictly  for  ihe  Electric 
fields  are  in  bold  face  and  commands  lhaf  are  strictly  to  accommodate  MACSYMA  are 
underlined . 

These  first  commands  are  in  the  preparation  stage  of  the  process. 

realpart(ex); 

imagpart(ex); 

subst(((((x-xl)*2>+((y-yi)‘2)+((z-zi)‘2))*(l/2)),r,R); 

Level  one: 

substCAr.(x-xi).8)$ 
substCA2\(y-yi),R)$ 
subst("A3\(z-zi),X)$ 


[#] 


1 


l  tv*  I  twt: 


£ 


£ 


i 

V, 

V. 


V, 


K 


substCB  r.“Ar‘2,x)i 
substCB2YA2"2,X)$ 
substCB3YA3"2,X)$ 

3Ub5tCB4-»-A2'.  txri*'A  1  *  *'A2'.X)t 
substCB4”"'A3‘,jxri  ""AC*  ”A3\X)$ 
substCB5”"”Ar.  ixri*'A2  "~AC.X)1 
substCBS”"”A3”.jxr!*”A2”»”A3”,X)$ 
substCB6'"”Ar,  lxr)*'A3  ‘»'AI'.X)t 
substCB6"*'A2”,jxrj*"A3"""A2‘,X)$ 
substCB7',"A2".  ixifAl”  »'A2'.X)t 
substCB7”"”A3\jxii*”Ar"'A3”,S)$ 
subst('B8‘" 'A  1”,  lxii**A2*  "'ArXli 
substCB8""”A3”,jxii"'A2'*'A3”,X)i 
SubStCB9""'Ar.  |xii*~A3~  »'AC.X)i 
substCB9‘ *  ”  A2” ,jxii*  "A3" » ' A2’ ,S)$ 
substCBlO”*'A2',  lvri"'A1  '"’A2VX)i 
substCB  10"*  "A3', jyri*  ”A  I  '•  "A3”,X)$ 
substCBI  r*"Ar,  ivri"'A2'»'AI~.X)l 
substCB  1  r*  ”A3\jyri*  ”A2”"  “A  3“.J8)$ 
substCB  12‘"  "AC,  |yrt""A3‘  #'A I '.%)! 
substCB  1 2"*  ”A2”,jyri*  'A3'*  -A2".J?)t 
substCBI 3'* ”A2”.  lvii"*A  C  »'A2'.^)t 
substCB  13'- 'A3'.jyUB'AC»'A3‘ 
substCB  1 4"  "‘AC,  ivii«‘A?-«-At-.X)l 
substCB  14'"  -A3".jyii  *  ”A2’"  *A3\X)t 
substCB  15"* ‘AC.  ivli"'A3*  »'A C  S)1 
substCB  15'-  "A2“.jyii  •  "A3""  ”A2\X)$ 
substCB 16'"'A2',  lzrl""Al~  »'A?'  X)1 
substCB  16'"  'A3”,jzri"”A  1  ”•  'A3'.X)$ 
substCB  17'* "AC,  izrt"*A2‘  ■*Ar.»i 
substCB  1 7"»  A3",jzri •  'A2“* "A3',S)t 
substCB  18'* 'AC,  lzrl*"A3”  *'AC.X)i 
substCB  18'""A2',jzn*  “A3”"  ”A2\X)$ 
substCB  19‘*”A2”,  lzii*”A  1  ‘  * "A?".X)t 
substCB  19”*  ”A3”.jzii  •  *A  C"  ”A3".X)| 
substCB20”*  AC.  izii»~A2~  *'AC.X)$ 
subst('B20“" 'A3”.jzu" 'A2”"  A3”.»)$ 
substCB2C"'AI",  izii»”A3”  " ‘A !”.%)> 
substCB2 1  ”•  ”A2”,jzii*  "A3”*  'A2”,f.)$ 

suk9t(*B22*.3ab*Jxri,%)$ 
subst(*B23*a3*b*  jxii,X)| 
subst(*B24\3»b»  jyrl,X)| 
S4ibstCB25\3*b*  jyll.X)$ 
siibst(*B26’,3*b*  jzrl,%)| 
subst(”B27‘.3"b"  JzllaX)t 


:%• 


subst('F9 \*E2*-  'B24*.X)t 
substCF  10*. *E2**’B25',JI)$ 
substCF  11 VE4-*  jyM.X)* 
subsICF  12*,’E4"“  jyll,X)| 


substCF  1 3 \*E2'»  *B26\X)t 
substCF  1 4VE2'-  *B27\%)$ 
substCF  l5VE4*«jzri.X)* 
substCF  16VE4’*  jzii,X)$ 


a i 


subst('620VF8'-*F5'.X)$ 
subst(*G2 1*,"F7"+*F6*,X)$ 
subSt('622','F4"*'D34',X)$ 
subs  t(  "623 VF4"  *  ‘035’, X); 

substC*630‘,‘F  1 2'-'F9*,X)$ 
subs  t( ’63 1  \'F 1 1  ■+*F10',X)$ 
substC“632*,'F4‘*  'D36\X)* 
subst(*633 \*F4'«  'D37*,X); 


subst('640‘.*FI6*-'FI3‘.X)$ 
subst("641‘,"F15’+"F14*,X)$ 
subsU'642 VF4"«  'D38\%)| 
subst(*643*,*F4‘*  “D39*,X); 


Level  eight: 


subst('H22  VC33'*  '62 1  \X)$ 
subst('H23*.'C33'«  '620'.*)$ 


sub3t('H32'.*C32'*  '63 1  '.%)$ 
subst(*H33’.'C32'«  *630'.%)$ 


Level  nine: 

subst(*J20\'H23'+'H2 1  *.*)$ 
subst(’J2 1  *.‘H22'*'H20\%); 
sgbst(*J30 \'H33'+'H3 1  '.*)* 
subst(*J3 1  VH32*+*H30*,%); 
sgbst('«i40'.*H43*+*H4 1  *.X)$ 
subst(*J41*>*H42'+*H40*,X); 


»l  »i  : »•  , 


Li  tkt  i-' 


4  (tl  *k*"  »v*  «  »  t.»«i,  «  >-.Mi  ^t.  ^*4 »■  .n  , 


.I*  ji  i.  ■.' _>.' 


Level  ten: 

substCK20  VF I  *«  *J2 1  *.X)$ 
subst(*K2 1  *,*F2‘«  *J20*,X); 
substCK22\‘F  I  **  -J20M)$ 
substCK23  VF2*«  M2 1  *.X); 

sub5t(’K30‘.*F  1*  •  *J3 1  \X)t 
subst(*K31","F2**“J30*,X); 
substCK32*.-F  t  "*  M30*.X)$ 
subst(-K33*.*F2'*  -J3 1  \X); 

subst('K40','F  1  *•  M4 1  *.X)$ 
subst(~K4 1  VF2*«  *J40‘.X); 
s«bstCK42*.*F  I '«  M40‘.X)$ 
substCK43  VF2**  M4 1  ‘.X); 


Level  eleven: 

subst(‘L20*,~K2 1  *+-K20\X); 
subst(*L2 1  *,"K22"-"K23“.X); 

substCL30  VK3 1  -+*K30*.X); 
subst('L3 1  VK32*-*K33\X); 

subst('L40“.’K4 1  ’♦-K40-.X); 
substCL41  '.'K42’-“K43’.X); 


Level  twelve: 


substCM  I  \*F3’«  *L20\X) 
subst(*t12*  ,'F3_*  *L2 1  *  .X) 
subst(’H3*,*F3"*  *L30'.X) 
subst(*M4  VF3‘*  *L3 1  ’.X), 
substCn5  VF3**  'L40\X); 
substCM6VF3*»*L4r.X) 


Wi 


This  appendix  provides  the  supplementary  material  for  the  develop¬ 
ment  of  the  Parallel  VWE.  Presented  here  are  the  full  size  data-flow  graphs 
which  provide  a  graphical  representation  of  the  specification  data  developed 
as  a  result  of  this  thesis  effort.  The  text  is  organized  into  three  sections. 

The  sections  are  the  presentation  of  the  Parallel  Vector  Potential,  the  Parallel 
Magnetic  Field,  and  the  Parallel  Electric  Field  The  order  of  presentation  for 
all  sections  are 

1.  The  real  part  of  the  "X"  directed  vector  quantity 

2  The  imaginary  part  of  the  vector  quantity 

3  The  real  and  imaginary  part  (paired)  of  the  vector  quantity 

4  The  combined  cartesian  set  for  the  vector  quantity 

The  final  data-flow  graph,  Figure  28,  is  the  presentation  of  the  fully 


overlayed  vector  quantities 
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